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On the eigenvalues for slowly varying perturbations of a periodic 

Schrodinger operator 

Magali Marx 


Abstract 

In this paper, I consider one-dimensional periodic Schrodinger operators perturbed by a slowly 
decaying potential. In the adiabatic limit, I give an asymptotic expansion of the eigenvalues in the 
gaps of the periodic operator. When one slides the perturbation along the periodic potential, these 
eigenvalues oscillate. I compute the exponentially small amplitude of the oscillations. 
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1 Introduction 


In this paper, we study the spectrum of one-dimensional perturbed periodic Schrodinger operators. Pre¬ 
cisely, we consider the Schrodinger operator defined on L^(R) by: 

Hip,e = -^ + [V{x) + W{ex + ip)], ( 1 . 1 ) 

where e > 0 is a small positive parameter, is a real parameter, and V is a real valued 1-periodic 
function. We also assume that V is and that W is a fast-decaying function. 

The operator can be regarded as an adiabatic perturbation of the periodic operator Hq-. 


Ho=-A + V. (1.2) 

The spectrum of the periodic operator Hq is absolutely continuous and consists of intervals of the real 
axis called the spectral bands, separated by the gaps. 

If the perturbation W is relatively compact with respect to Hq, there are in the gaps of Hq some 
eigenvalues EHIEI] We intend to locate these eigenvalues, called impurity levels. 

The equation 

( 1 - 3 ) 


depends on two parameters e et ip. We study the operator in the adiabatic limit, i.e as e ^ 0. The 
periodicity of V implies that the eigenvalues of are e-periodic in ip. We shall shift ip in the complex 
plane and we shall assume that W is analytic in a strip of the complex plane. 

If V = 0, there are many results. The case when W is a well has been studied; in the interval ] inf W, 0[, 

M. 

there is a quantified sequence of eigenvalues [3] . We shall give an analogous description of the eigenvalues 
of in an interval J out of the spectrum of Hq. Precisely, when W and J satisfy some additional 
conditions described in sections r2.2.1l I2.2!31 et 12.41 we show that the eigenvalues of oscillate around 
some quantized energies. The quantization is given by a Bohr-Sommerfeld quantization rule; the ampli¬ 
tude of oscillation is exponentially small and is determined by a tunneling coefficient. 
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1.1 Physical motivation 

The operator is an important model of solid state physics. The function 'ip is the wave function 
of an electron in a crystal with impurities. V represents the potential of the perfect crystal; as such 
it is periodic. The potential W is the perturbation created by impurities. In the semiconductors, this 
perturbation is slow-varying m- It is natural to consider the semi-classical limit. 

1.2 Perturbation of periodic operators 

In the spectral theory of the perturbations of a periodic operator 


Hp = Ho + P (1.4) 

has motivated numerous studies with different view points. 

The characterization of the existence of eigenvalues is not easy: particularly, in any dimension, m deals 
with the existence of embedded eigenvalues in the bands. On the real axis, the situation is simpler. When 
the perturbation is integrable, the eigenvalues are necessarily in the adherence of the gaps (Ealia). 

To count the eigenvalues in the gaps, many results have been obtained thanks to trace formulas. In the 
large coupling constant limit, i.e when P = XU, with A — *■ -|-oo, I2S1IIIE1 have studied lim tr(p/^^^,,), 

where = ^[e,e']H\ (spectral projector of H\ on an interval [E,E'] of a gap of PIq). In the semi- 

classical case, ^ has given, under assumptions close to mine, an asymptotic expansion of tr[/(iJ,p_£)], for 
/ S (^“(IR.) and Supp / in a gap of Hq. These formulas are valid in any dimension but are less accurate. 
For example, in the expansion obtained in the accuracy depends on the successive derivatives of the 
function /; the formula does not give an exponentially precise localization of the eigenvalues. 

In the one-dimensional case, the scattering theory, well-known in the case V = 0, has been developed in 
cniiiTi for the periodic case. Precisely, we construct some particular solutions of equation (iia, which 
tend to zero as x tends to infinity. We call these functions recessive functions. The eigenvalues of equation 
are given by a relation of linear dependence between these solutions. 

1.3 Main steps of the study 

We give here the main ideas of the paper. An important difficulty is the dependence of the equation on 
the parameters e and v?; particularly, one has to decouple the “fast” variable x and the “slow” variable 
ex. The new idea developed in PEI is the following : we construct some particular solutions of 
satisfying an additional relation called the consistency condition: 

f{x + l,ip,E,e) = f{x,ip + e,E,e). (1.5) 

This condition relates their behavior in x and their behavior in ip. 

To hnd a recessive solution of it suffices to construct a solution of which satisfies and 

which tends to 0 as |Re ip\ tends to -|-oo. First, we build on the horizontal half-strip € C ; ip & 

] — oo, —A\ -I- i[—y, y]} a solution h^_ of equation HI ..‘111 which is consistent and which tends to 0 as Re 
tends to — oo. Similarly, we construct for {(^ G C ; p € [A, -|-oo[-|-f [—y, y]} (Theorem |2l. These 
functions are recessive for the variable x. The characterization of the eigenvalues is given by the relation 
of linear dependence between h® and : 

w{h^_, = 0. 

In the above-mentioned equation, w represents the Wronskian whose definition is recalled in (lO) . 

It remains to compute w{h^_, h’P^). To do that, we use the complex WKB method developed by A. Fedotov 
and F. Klopp. This method consists in describing some complex domains, called canonical domains, on 
which we construct some functions satisfying dni) and having a particular asymptotic behavior: 

f±{x,p,E,e) = '^{'ip±{x,p,E) +o{l)), e ^ 0. (1.6) 
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In equation Oil . the function k is a analytic multi-valued function, defined in TIM : the functions ili± 
are some particular solutions of equation 


Ho^ = {E-Wi:p))^, 

analytic in (p on these canonical domains and called Bloch solutions. We will prove the existence of such 
functions in section E21 

A. Fedotov and F. Klopp prove the existence of functions with standard asymptotic only on compact 
domains of the complex plane. We shall extend some results on infinite strips of the complex plane. The 
consistency condition implies that the function h® satisfies the standard asymptotic HI .fi|l to the left of 
—A and that h'^ satisfies an analogous property to the right of A. Thus, the computation of w{h^_, h‘^) 
is similar to the calculations of A. Fedotov and F. Klopp. We must find a sufficiently large domain of the 
complex plane, in which we know the Wronskian of h® and h'^. 

The methods used in their works underline some topological obstacles, which change the standard asymp¬ 
totic (II. til) : these obstacles depend on W and E. We give precise assumptions in sections TIM and TIM 

2 The main results 

In this section, we describe the general context and the main results of the paper. 

First, we present the assumptions on the potentials V and W, and on the interval J. There are mainly 
three kinds of assumptions. Firstly, the study requires some assumptions on the decay of W to develop 
the scattering theory. Then, in view of the hypotheses of the complex WKB method of [HI, we assume 
that W is analytic in some domain of the complex plane. Finally, we shall depict the geometric framework 
and particularly the subset [E — 

We obtain an equation for the eigenvalues in terms of geometric objects depending on i?o, ^ and E: the 
phases and action integrals, defined in sections ESI 

2.1 The potential V 

We assume that V has the following properties: 

(Hv.p) V is 1-periodic. 

We consider as a perturbation of the periodic equation: 

+ V{x)iljix) = {E -W{p))tp{x). (2.1) 

We shall use some well known facts about periodic Schrodinger operators. They are described in detail 
in section 0 

We just recall elementary results on Hq. The operator Hq defined in is a self-adjoint operator on 
71^ (R). The spectrum of Hq consists of intervals of the real axis: 

O’iHo) = [^[E2n+l,E2n+2], ( 2 . 2 ) 

neN 


such that: 

El < E 2 < Eq < Ei...E2n h E2n+1 < E2n+2-", En —*■ +00,71 —*■ -fOO. 

These intervals [if 2 n-i-i, E 2 n+ 2 ] are called the spectral bands. We set Eq = —oo. The intervals (i? 2 „, E 2 n+i) 
are called the spectral gaps. If E 2 n ^ E 2 n+i, we say that the gap is open. 

Furthermore, we assume that V satisfies: 

(Hv,g) Every gap of Hq is not empty. 
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This assumption is “generic”, we refer to m section XIIL16. An important object of the theory of 
one-dimensional periodic operators is the Bloch quasi-momentum k (see section El. This function is a 
multi-valued analytic function; its branch points are the ends of the spectrum, they are of square root 
type. We shall give a few details about this function in section^ Finally, we suppose: 

(Hv) V satisfy (Hv,p) and (Hv,g). 


2.2 The perturbation W 

2.2.1 Smoothness assumptions 

We assume that W is such that: 


(Hw.r) There exists Y > 0 such that W is analytic in the strip Sy = {|Im (^)| < Y} and 
there exists s > 1 et C > 0 such that for z G Sy, we have: 


|W(z)| < 


C 

1 + |zh’ 


(2.3) 


These assumptions are essential to develop the complex WKB method. The analyticity of the perturbation 
is crucial in the theory of [S] . The decay of W replaces the compactness resulting from periodicity in El ■ 
We begin with presenting the complex momentum. This main object of the complex WKB method shows 
the importance of W“^(R). 


2.2.2 The complex momentum and its branch points 

We put: 

C+ = {ip G C ; Im > 0} and C_ = {yj S C ; Im ip < 0}. 

For equation Oil . we consider the analytic function k defined by 

k(p) = k(E — W(p)}. (2.4) 

We recall that the function k is presented in section im The function n is called the complex momentum. 
It plays a crucial role in adiabatically perturbed problems, see ElEl. 

N is the set of non-negative integers. We define: 

T(£;) ={pgSy ■, 3nGW / E- W{p) = £!„} (2.5) 

The set of branch points of k is clearly a subset of T {E). The following result gives a characterization of 
the branch points of n among the points of T {E ): 

Lemma 2.1. Let p be a point of T{E). //inf{g ; 0} G 2N + 1, then p is a branch point of k. 

This result follows from the fact that the ends of the spectrum are of square root type. 


2.2.3 Geometric assumptions 

The spectrum a{Ho) consists of real intervals. Fix E G M.. If E — W{p) is in the spectrum a{Ho), then 
W(p) is real. The spectral study of 111.311 is then tightly connected with the geometry of W“^(R). 

We state now the geometric assumptions for W. These assumptions are mainly a description of VF“^(R) 
in a strip containing the real axis. We call strictly vertical a line whose slope does not vanish; for precise 
definitions, we refer to section Em 
(Hw.g) 

1. W|R is real and has a finite number of extrema, which are non-degenerate. 
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Figure 1: A subset of W ^(R) 


2. There exists Y > 0 and a finite sequence of strictly vertical lines containing an ex¬ 
tremum of W, such that: 

W-^(M)nSY= U (Si)UM. (2.6) 

ie{i...p} 

2.3 Some remarks 

• Since W is real analytic, we know that W{jf) = W{}p)\ this implies that VF“^(R) is symmetric with 
respect to the real axis. 

• We define =2^(1 C+ and = 11^(1 C_. 

• Figure n shows an example of the pre-image of the real axis by such a potential. 

As we have explained in section o we cover the strip Sy with local canonical domains. On these 
domains, we construct consistent functions with standard behavior (ie satisfying dca and ( 11 . 61 ) 1 . 

To compute the connection between the bases associated with different domains, we get round the 
branch points (for analog studies, we refer the reader to ^2 El)- We will now state some more ac¬ 
curate assumptions about the configuration of the branch points; in particular, these assumptions specify 
{E — W)~^ {a{H[))) when E is real. The spectral results of A. Fedotov and F. Klopp on perturbed periodic 
equation have shown the importance of the relative positions of J and cr(iJo)- 


2.4 Assumptions on the interval J 

Now, we describe the interval J on which we study equation in. 

2.4.1 Hypotheses 

We assume that the interval J is a compact interval satisfying: 

(Hj) 
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Figure 2: [E—W) 


1. For any E S J, there exists only one band B of o-(Ho) such that the pre-image C := (E — W) ^(B) 
is not empty. 

O 

2. For any E € J, C := (E ~ W)~^(B) is connected and compact and (E — W)“^(b) contains 
exactly one real extremum of W. 

2.4.2 Consequences 

• [Hj) implies that J is included in a gap. 

• The band B in [Hj) (1) depends a priori on E. But, since J is connected, the band B is fixed for 
any E G J. 

• Similarly, the extremum of W in assumption (Hj) (2) depends on E, but by connectedness, it is 
the same for any E £ J. 

2.4.3 Notations 

Put B = [E 2 n-i, E 2 n], for n G N*. Moreover, we can always change VF or so that the extremum of W 

in (2) is 0. 

Then {Hj) has the following consequences: 

1. For any E G J, {E - W)-\a{Ho)) n Sy = {E - W)-\B) n Sy 

2. Let Ej. G {E 2 n-i, E 2 n} be the end of B satisfying E^ G {E — 1F)(]R) for any E G J. We define E^ 
such that {E^,Er} = {E 2 n-l,E 2 n}- 

3. There are exactly four branch points {(fir , ‘fit ) ^ fo related to Er and Ei. They 

satisfy: 

E-W{ipt) = Er,E-Wi^-) = Er, ip- <0<(^+ 

E - W{pi) = E - W{^) = Ei, Im Pi > 0. 

4. There exists a strictly vertical line a containing 0 and connecting pitopi, such that {E—W)-^{B)r] 

Sy = [p- ,Pr\^<J. We define cr+ = (TnC+ and ct_ = aOC-. We let E = (E—IF)“^(R)\R, ct C E. 

These objects are described in figure |2 
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Figure 3: Some examples of potential W 


2.4.4 Remarks and examples 

We first give a few comments on assumption (ffj). 
• We call C the cross. 


• This assumption means intuitively that, in Sy, we see the band B only near the extremum 0. 

To illustrate these technical assumptions, we give a few examples of potentials W and intervals J. We 
have depicted some examples in figure |21 

• The simplest case is when W has only a non-degenerate minimum W- (see figure 01 A), 
in concrete terms, we can think of the example: 

2 ’ a>0, 

1 + 

Then, if we fix i? = [E 2 n-i,E 2 n] and T < 1, we can choose J = [a, b] such that: 

OC Q 

max{F; 2 n- 2 , E 2 n -1 - a, E 2 n - ^_y 2 ^ < a <b< min{E 2 n-l,E 2 n - a, E 2 n +1 - 


• We can assume that W has a maximum W+ and a minimum W_, if J is chosen to see the band 
only near the maximum (see figure OjB). 


W{x) 


2 1 
1 + x^ ~ 1-f (a; - 5)2 


J C]i?2n-1 + i?2„_2 + [U]i?2„, £^271-1-1 + [, |J| < |£^271-2 — £- 2 , 1 - 1 1 


Consider this example a little further. The choice of Y is more complicated in this case. The 

study of equation W{u) = w for w > W+ shows that there exists only one solution in the strip 

{Im u €]0, ![} that we call Z(w) ; we choose V €] sup Z(E — E 21 - 1 ), inf Z(E — E 2 i- 2 )l- 

EeJ EeJ 


• In fact, we could adapt our method to weaker assumptions. For example, we can assume that we 
do not see the branch points ipi and (incomplete cross), which means that the vertical line cr 
does not contain any branch points of k. We refer to section 1^.5.41 for some details. 


• For the sake of simplicity, we have assumed that all the extrema of W are non degenerate. Actually, 
it suffices to assume that only the extremum of W in 0 is non degenerate. 

• Similarly, we could weaken assumption {Hv,g)- We only have to assume that the gaps adjoining 
the band B of (i£j) are not empty. 
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2.5 Phases and action 


In this section, we define the tunneling coefiicient t and the phases $ et these analytic objects play 
an essential role in the location of the eigenvalues. These coefficients are represented as integrals of the 
complex momentum k in the ip plane. 

In the strip Sy, we consider k a branch of the complex momentum, continuous on C. 


2.5.1 Definition and properties 

We introduce the action S and the phases <1> and related to the branch k. 
Definition 2.1. We define the phase: 


the action: 


$(£;) = 

fvt 

/ {k(u) - K{ipf))du, 

(2.7) 


'vr 


S{E) = 

i / ( k ( m ) - K{ifi))du, 

J a 

(2.8) 


the second phase: 


= / {K{u) — K{if^))du+ {K{u) — K{(pf))du+ {K{u) — K{(pi))du— {k{u) — K{ipi))du. (2.9) 
d tfP J {ff J d 

In section IHl we prove the following result on the behavior of the coefficients $, S and d)^. 

Lemma 2.2. There exists a branch ki such that the phases and action integrals have the following 
properties: 


1. S, o,re analytic in E in a complex neighborhood of the interval J. 

2. S, take real values on J. and S are positive on J. 

3. WE € J, ^'{E){E, -Er)>0, SiE) < 2 tt Im {ipi{E)). 

We define the tunneling coefficient : 

t{E,e) = eyipi{—S{E)/e). (2.10) 


t is exponentially small. 

2.5.2 Remark 

The phase and action are simply a generalization of the coefficients of the form / ^JE — W{ip)dip, well- 
known in the case V = Q (we refer to ismniiHi) 

We point out that the coefficient <!> depend only on the value of W on the real axis, whereas S and 
depend on the values of W in the complex plane. The phase is independent of the analyticity of W 
unlike S and 

Now, we state the equation for eigenvalues for 1131). 

2.5.3 The main result 
Theorem 1. Equation for eigenvalues. 

Let V, W and J satisfy assumptions {Hy), (H\y,r), [Hw.g) OLnd [Hj). FixY^ g]0, T). 

There exists a complex neighborhood V of J, a real number Eq > 0 and two functions 3> and ^d with 
complex values, defined on Vx]0,eo[ such that: 



• The functions ‘&(-,e) and are analytic on V. Moreover, and satisfy: 

^{E,e) = ^{E) + ho{E,e) and ^d{E,e) = ^d{E) + hi{E,e), 
where p is a real coeffieient, ho{E,e) = o(e) and hi{E,e) = o(e) uniformly in E & V. 

• If we define the energy levels in J by: 

^EW{e),e) V/€{L_(e),...,L+(e)}, (2.11) 

then, for any e G]0,eo[, 

— the spectrum of in J consists in a finite number of eigenvalues, that is to say 

U {Ei{p,e)}, (2.12) 

iG{Z/_(e)....,L+(e)} 


— these eigenvalues satisfy 


where there exists c > 0 such that 


( (e), e) + 27 r(j 3 + pe 


(2.13) 


sup r{E,(p,e) < —e 
£;gv.¥>gr c 


We prove this result in section |S1 


2.5.4 Remark 

If we only assume that cr does not contain any branch points, asymptotic is replaced by the 

estimate: 

\Ei{ip,e)-E^^\e)\<Ce-^ 
where 2Y is the width of the strip Sy- 

2.5.5 Application : asymptotic expansion of the trace 

By using the previous result, we can compute the first terms in the asymptotic expansion of the trace 
formula , and partially recover a result of 0]. 

Corollary 1. Let f G (^“(R) be a real function such that Supp f G J. Then the function f{Hp,^g) is 
e-periodic in p and its Fourier expansion satisfies: 

tr [f{Hp,g)] = - nf{Hu,e)]du + 0(e-®/") (2.14) 

£ Jo 

[ [f {Hu,e)]du =-^ [ [ f{W{u) + E{K))dKdu + o{e) (2-15) 

Jo ^ 7 ’" Jr„ JI-tv,-k] 

where S = inf ^(e) >0 

eGSupp f 

We give more details and the proof of this corollary in section I8.fil 


9 







3 Main steps of the study 

Here, we explain the main ideas of the paper. 


3.1 One-dimensional perturbed periodic operators 

3.1.1 

We consider equation dOl) as a perturbation of the periodic equation 

Ho^J = (3.1) 

where the operator Hq is defined in To do that, we shall describe the spectral theory of periodic 

operators in sectional 

• For the moment, we simply introduce the Bloch solutions of equation EU. We call a Bloch solution 
of (12.111 a function dt satisfying (12.111 and: 

Va;GK, '^ix + l,E) = X{E)'^{x,E), (3.2) 

with A 0 independent of x. The coefficient X{E) is called Floquet multiplier. We represent X{E) in 
the form X{E) = k is the quasi-momentum presented in section |0] and described in section 

0] If if ^ (T(ilo), there exist two linearly independent Bloch solutions of Il3.1ll fsee section ICTl . We 
call them et dt-; the associated Floquet multipliers are inverse of each other and the functions 
dtj- are represented in the form: 

'^±[x,E) = e^'^^^^'^^p±[x,E) avec p±{x+ 1,E) = p±{x,E). 

For Im k{E) > 0, the function d'-|_(a;,if) tends to 0 as a; tends to -l-oo and the function ^-{x,E) 
tends to 0 as a; tends to — oo. Actually, equation (|S3 defines the functions d'+ and dt- except for 
a multiplicative coefficient. Precisely, equation E3 defines two one-dimensional vector spaces that 
we call Bloch sub-spaces. 

To study the eigenvalues of perturbations of periodic operators, nni and na introduce, for Im k{E) > 
0, two functions {x,ip,E,e) i—> F+{x,ip, E,e) and {x,ip,E,e) F-{x,(p, E,s) solutions of lESl 
satisfying: 

lim [F+{x,(p,E,e) — 'i'+{x,E)] = 0, lim [F-{x,ip, E,e) —'i’-{x, E)] = 0 (3-3) 

Condition (Tnil guarantees the uniqueness of F+ (resp. of FT) since the function d'+ (resp. d^-) 
tends to 0 as a: tends to -l-oo (resp. —oo). These functions are called dost functions; they are 
generally constructed as solutions of a Lippman-Schwinger integral equation. This construction is 
an adaptation of the usual theory of scattering (chapter XI of m) for a perturbation of laplacian; 
it consists in looking for particular solutions of ESI from the solutions of the periodic equation. 

We call Jost sub-spaces the sub-spaces f7+ and f/b generated by F+ and FA. 

J+ (resp J-) is the set of solutions of HI.311 being a member of F^([0, oo)) (resp. F^((—oo, 0])). 

• Let / and g be two derivable functions, the Wronskian of / and g called w{f,g) is defined by: 

w{f,9) = f9-f9 (3-4) 

We recall that if / and g are the solutions of a second-order differential equation, their Wronskian is 
independent of x. The spectral interest of the Jost sub-spaces is the following: 

Proposition 1. We assume that Im k{E) > 0. Let h~ G and h'^ G J7+ be two nontrivial Jost 

solutions of E3. E is an eigenvalue of if and only if: 

w{h+,h~)=0 (3.5) 

To compute the eigenvalues, it suffices to construct the Jost sub-spaces. 
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3.2 Construction of consistent Jost solutions 

We denote by (Hj) the following assumption: 

(H5) There exists n G N such that J is a compact interval of ]E 2 n,E 2 n+i[- 

Clearly, {Hj) is weaker than (Hj). 

We introduce a new notation. 

For a function f : U C C” ^ C^, we define the function /* : U ^ 

nZ)=W)- (3.6) 

As we have explained in section ^ an useful idea to study dOJ is the construction of consistent solutions, 
i.e. satisfying dinii. First, we choose in and J7+ some consistent bases. We shall prove the following 
result: 

Theorem 2. We assume that {Hy), {Hw,r) o-nd (Hj) are satisfied. Fix X > 1. Then, there exist a 
complex neighborhood V = V of J, a real Eq > 0, two points mg and md in C, two real numbers Ag and 
Ad and two functions {x,(p,E,e) h^_{x,(p, E,e), {x,ip,E,e) i—!■ hi^{x, (f, E, e) such that: 

• The functions {x,(p,E,£) h^{x, ip, E, e) and {x,ip,E,£) i—> h'^{x,ip,E,e) are defined and consistent 

onRx Sy X Vx]0, eo[. 

• For any X G [—X, X] and£G]Q,£o[, {(fi, E) h^_{x, ip, E, e) and {p, E) h'^{x,p, E,£) are analytic 
on Sy X V. 

• The function X h^{x, p, E, e) (resp. x i—> h'^(x, p, E, e)) is a basis of ff- (resp. J+). 

• The functions hf and h‘^ have the following asymptotic behavior: 

h^- (x, p, E, e) = e~ ■^’” 9 .0_ (^x, p, E){1 + Rg{x, p, E, e)), (3.7) 

and 

h'lix, p, E, e) = e^ p, E){1 + Rd{x, p, E, e)), (3.8) 

where 

— Rg and Rd satisfy: 


sup max{\Rg{x,p,E,E)\,\da;Rgix,p,E,E)\} <r{E). 

xE]-X,X[, Re v<Ag,E^V 

sup m.ax{\Rd(x,p,E,E)\, \d,:Rd(x,p,E,E)\} < r(e:), 

xE]-X,X[, Re ip>Ad,EEV 

with 

lim r{E) = 0. 

— The functions tp^ and ip- are the Bloch canonical solutions of the periodic equation eu defined 
in section \5.A 

• There exist two real numbers Og G { — 1, 1}, Od G { — 1, 1}, an integer p and two functions E 1 —s- ag(E) 
and E 1 —!■ ad(E) such that: 

1. For any e G]0,eo[) x G M, E G V,et p G Sy ,we have: 

a*giE)(h<L)*{x,p,E,£) = iage-i^P^^ag(E)h^_ix,p,E,E) (3.9) 


ad(E)(h‘\f)*{x,p,E,£) =icTde'^^P^'^ad(E)h'f^(x,p,E,£) (3.10) 

2. The functions ag and ad are analytic and given by and (16.251) . They do not vanish on V. 
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We immediately deduce from Theorem |21 and Proposition Q] that the eigenvalues of are char¬ 
acterized by: 

= 0 (3.11) 

Theorem 12 is the consequence of two main ideas: 

• First, we adapt the construction of dost functions developed by nnmni. Indeed, this construction 
proves that asymptotic IH.lll is only satisfied on domains which depend on e (see section EJ- 

• We must understand how this asymptotic evolves on a domain which does not depend on e. To do 
that, we extend the continuation results of Fedotov and Klopp in a non compact frame (see section 

0 . 


3.3 Conclusion 

To finish the computations, it suffices to apply the methods of [S]- We have to go through the cross (see 
figure 0. We will show that there exists in the neighborhood of the cross a consistent basis with 
standard asymptotic. We shall express the functions /i® and h'^ on this basis (section 0. 


4 Periodic Schrddinger operators on the real line 

We now discuss the periodic operator where F is a 1-periodic, real-valued, Lf^^-function. We collect 
known results needed in the present paper (see nnnHnni). 

4.1 Bloch solutions 


Let 41 be a solution of the equation 

iJo^' = £4' 

(4.1) 

satisfying the relation 

4i(a: + 1) = A4'(a;), VgK 

(4.2) 


for some complex number A 7 ^ 0 independent of x. Such a solution is called a Bloch solution, and the 
number A is called the Floquet multiplier. Let us discuss the analytic properties of Bloch solutions. 

In 12.211 . we have denoted by [£li, £^ 2 ],..., [£l 2 n-i-i, £’ 271 + 2 ], • ■ • the spectral bands of the periodic Schrddinger 
equation. Consider r± two copies of the complex plane £ G C cut along the spectral bands. Paste them 
together to get a Riemann surface with square root branch points. We denote this Riemann surface by 
F. 

One can construct a Bloch solution 4(x, S) of equation 14.111 meromorphic on F. The poles of this solution 
are located in the spectral gaps. Precisely, each spectral gap contains precisely one simple pole. This 
pole is situated either on r+ or on r_. The position of the pole is independent of x. For the details, we 
refer to 0 J]. 

Except at the edges of the spectrum (i.e. the branch points of F), the restrictions 4+ of T on r± are 
linearly independent solutions of l4.ll Along the gaps, these functions are real and satisfy : 

4±(a;,£ - lO) = $±(a:,£-I- to), V£ €]£ 271 , £ 277+1 G N. (4.3) 

Along the bands, we have : 

T±(a:,£ - iO) = $=p(x,£-I-iO), V£ g]£27i+i, £ 277 + 2 ], « G N. (4.4) 
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4.2 Bloch quasi-momentum 

4.2.1 

Consider the Bloch solution 4'(a;,£) introduced in the previous subsection. The corresponding Floquet 
multiplier X{£) is analytic on F. Represent it in the form: 

X{£) = exp{ik{£)). (4-5) 

The function k{£) is called Bloch quasi-momentum. It has the same branch points as \E'(a;,£), but the 
corresponding Riemann surface is more complicated. 

To describe the main properties of k, consider the complex plane cut along the real line from Ei to -l-oo. 
Denote the cut plane by Cq. One can fix there a single valued branch of the quasi-momentum by the 
condition 

iko{£) <0, £ < El. (4.6) 

All the other branches of the quasi-momentum have the form ±/co(^) + 27rm, m € Z. The ± and the 
number m are indexing these branches. The image of Cq by ko is located in the upper half of the complex 
plane. 

Im koi£) >0, £g Co. (4.7) 

In figure 21 we drew several curves in Cq and their images under transformation E i—> ko{E). The quasi¬ 
momentum ko{E) is real along the spectral zones, and, along the spectral gaps, its real part is constant; 
in particular, we have 

ko{Ei)=0 ko{E2i±iO) = ko{E2i+i±iO)=±nl, IgN. (4.8) 

All the branch points of k are of square root type. Let Em be one of the branch points of k. Then, each 
function: 

/±(£) = {ko{£ ± iO) - koiEm ± iO))/V£ - Em, E gR (4.9) 

can be analytically continued in a small vicinity of the branch point Em- 
Finally, we note that 

ko{£)=V£ + 0{l/V£), l^l^oo (4.10) 

where E G Cq and 0 < aigE < 2 tt. 

The values of the quasi-momentum fcg on the two sides of the cut [Ai,-|-oo) are related to each 
other by the formula: 

y£G]Ei,+oo[, ko{£ + i0) = -ko{£-i0), Ei<£. (4.11) 

Consider the spectral gap (i? 2 /j A' 2 ;+i), I G N. Let C; be the complex plane cut from —oo to E 21 and from 
E 21+1 to -|-oo. Denote by ki the branch of the quasi-momentum analytic on C; and coinciding with fep 
for Im E > 0. Then, one has: 

y£ G] — 00 , E 2 i[U]E 2 i+i,+oo[, ki{£ + iO) + ki{£ — iO) = 2Trl. (4-12) 


4.3 Periodic components of the Bloch solution 

Let D a simply connected domain that does not contain any branch point of k. On D, we fix an analytic 
branch of k. Consider two copies of D, denoted by D±, corresponding to two sheets of Q. Now we redefine 
'I'± to be the restrictions of 4* to D±. They can be represented in the form: 

^±{x,£) =e^*’^^^^^p±{x,£), £gD (4.13) 

where p^{x,£) are 1-periodic in x, 

p±{x + 1,£) = p±{x,£), Va; G K (4-14) 
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Figure 4: The quasi-momentum k 


4.4 Analytic solutions of (j4.1|) 

To describe the asymptotic formulas of the complex WKB method for equation one needs specially 

normalized Bloch solutions of the equation (EH. 

Let 27 be a simply connected domain in the complex plane containing no branch point of the quasi¬ 
momentum k. We fix on 27 a continuous determination of k. We fix Sq G 27. We recall the following 
result (ISlini)- 

Lemma 4.1. We define the functions g± : 


and the functions : 


fg Pzp(x,£)dsp±(x,£)dx 
fgP+(x,£)p-(x,£)dx 


(4.15) 


-0^ : Kx27^C; (x,£) ^ (4.16) 

The functions £ 1 -^ '0^(x,f) are analytic on 27, for any a; G K. The functions '0^ are called analytic 
Bloch solutions normalized at the point £q of EIJ. 


Sometimes, we shall denote 0^(x,f,£o) to specify the normalization. We refer to section 1.4.4 of 
for the details of the proof. The proof follows from the study of the poles of 4^ and the zeros of k'. 
The poles of g± are simple and exactly situated at the singularities of v^4'±- The computation of the 
residues of g± at these points completes the proof. 


4.5 Useful formulas 

We end this section with some useful formulas. We recall that the functions g± are given in 14.1511 . Fix 
n G N. Equations and ESJ) lead to the following relations: 

gli£)=g±ix,£), y£ e]E 2 n , E2n+l [. (4.17) 

gl{£)=g^{x,£), y£ G\E2n+uE2n+2[. (4.18) 
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5 Main tools of the complex WKB method 

In this section, we recall the main tools of the complex WKB method on compact domains. The idea 
of the method is to construct some consistent functions of (II.dll with asymptotic behavior HI.(ill . This 
construction is not possible on any domain of the complex plane but on some domains called canonical. 
We apply the results of iniEip to the assumptions [Hw,g) and {Hj). We build a neighborhood of the 
cross, in which we construct a consistent basis with standard behavior In this section, we fix Y 

such that the assumptions {Hw,g) and {Hj) are satisfied in the strip Sy- 

5.1 Canonical domains 

The canonical domain is the main geometric notion of the complex WKB method. 

5.1.1 The complex momentum 

The canonical domains can be described in terms of the complex momentum Remind that this 

function is defined by formula iim . We have described k in section The properties of k depend 

on the spectral parameter E and of the analytic properties of W. 

We first formulate some definitions ([HI)- 

5.1.2 Vertical, strictly vertical curves 

Definition 5.1. We say that a curve 7 is vertical if it intersects the lines Im z = Const at non-zero 
angles 9. 

We say that a curve 7 is strictly vertical if there is a positive number 5 such that, at any point of 7 , the 
intersection angle 9 satisfies the inequality: 

S < 9 < TT — 6. (5-1) 


5.1.3 Canonical, strictly canonical curves 


Let 7 be a vertical curve which does not contain any branch point. On 7 , fix a continuous branch of the 
momentum of k. 


Definition 5.2. We call 7 canonical if, along 7 , 

• Im If I—> Im K{u)du is strictly increasing. 

• Im (p 1 -^ Im {k{u) — Tr)du is strictly decreasing. 

Assume that 7 is strictly vertical. If there is a positive number 5 such that, along 7 : 


rv' 

Im / K,{u)du > (5Im [p' 
J ip 


‘p) G 7^ 


(5.2) 


and 


j-v' 

Im / (tt — K{u))du > 5Im {p' 
J ip 


P) V(v?,^p')G7^ 


(5.3) 


we call 7 (5 — strictly canonical. 

We identify the complex numbers with vectors in To construct canonical lines, we have to study the 
vector fields n and k — tt, or rather their integral curves. For p £ D, S{p) denotes the sector of apex p 
such that, for any vector z G S{p), we have: 


Im {iHi(p){z — p)) > 0 et Im {i{K{p) — n){z — p)) < 0. 


(5.4) 


Let 7 S D a curve which does not contain any branch point. For all p G we denote t{p) the vector 
tangent to 7 in and oriented upward. The curve 7 S D is canonical for the determination k if and 
only if for any p G "f, the vector t[p) belongs to S{p) (see figure 0. The cone S{p) depends on the 
determination of k. For example, if k satishes Re k s]0, tt), this cone is not empty. 
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Figure 5: The cone S{ip) 


5.1.4 

In what follows, .fi and ^2 are two points in C such that 

Im < Im ^ 2 - 

We shall denote by 7 a smooth curve going from to ^ 2 ', this curve will always be oriented from to 

6 - 


5.1.5 Definition of the canonical domain 

Let iF be a simply connected domain in {Im ip G [Im ^i,Im ^ 2 ]} containing no branch points of the 
complex momentum. On K, fix a continuous branch k. 

Definition 5.3. We call K a canonical domain for k, and ^2 if it is the union of curves that are 
connecting and ^2 and that are canonical with respect to k. 

If there is 6 > 0 such that K is a union of S-strictly canonical curves, we call K 5 — strictly canonical. 

5.1.6 

Assume that AT is a canonical domain. Denote by dK its boundary. Fix a positive number 6. We call 
the domain 

C = {z € K ; dist{z,dK) > (5} 

an admissible sub-domain of K. 

Note that the branch points of the complex momentum are outside of C, at a distance greater than S. 

5.2 Canonical Bloch solntions 

To describe the asymptotic formulas of the complex WKB method for equation we shall use the 

analytic Bloch solutions of defined in Lemma EU for the parameter £ = E — W{(p). Precisely, we 

consider the unperturbed periodic equation: 

Hofj = {E-W{ip))f,. (5.5) 


5.3 

Let D be a simply connected domain in Sy, containing no branch points of k. The mapping ip 1 —> E—Wtpp) 
maps D onto a domain V G C. The domain V does not contain any branch point of k. 
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Fix ifQ € D, such that k'{E — W{‘f[j)) ^ 0. In Lemma f4.II we have built the analytic Bloch solutions 
{-0^} of equation (14.Ill , normalized hi E — VF((/5o)- For ip G D, we define: 

tjj±{x,ip, E) =E — W{ip)), Vm G R, Wip G D. (5.6) 

In 1^, it is proved that the functions p i—> ^p±(x, p, E) can be analytically continued to D. ^p± are called 
the canonical Bloch solutions of equation (1^ . Sometimes, we shall precise ip±{x,p, E,pq) to specify 
the normalization. 

We define 

u;±{p,E) = -W\p)g±iE-Wip)). (5.7) 

We also define: _ 

q(p) = ^k'(E - W(p) (5.8) 

5.4 The consistency relation 

5.4.1 Consistent functions and consistent bases 

We recall that we say that / is a consistent function if it satisfies deg. 

We say also that a basis {/±} of solutions of (II.dll is a consistent basis if: 

• The functions /+ and /_ are consistent. 

• Their Wronskian is independent of p. 

5.4.2 Analyticity and consistency 

First, we define the width of a set. 

Definition 5.4. Fix lb > 0 and M C Sva « set of points. We define 1 (M,Yq): 

1{M, Yq) = inf sup {\Re p — Re p'l, (p, p') G such that Im p = Im p' = y| (5.9) 

yG[-Yo.Yo] 

1(M,Yq) is called the width of M in Syq 
One has: 

Lemma 5.1. Fix E. We consider X > 0, Y G]0,y[, Eq > 0 and K a complex domain such that 
l(K,Y) > Sq. We assume that for any e g] 0, eo[; /(g i?, e) is a consistent solution of (II.dll for p G K 
and that for any x G [—X,X], the function p f(x,p,E,e) is analytic on K. Then, for any e G]0,eo[ 
and any x G [—X,X], the function p f(x,p,E,e) is analytic on Sp. 

This result is proved in 130. 

5.5 The theorem of the complex WKB method on a compact domain 

In this section, we recall the main result of the complex WKB method. 

5.5.1 Standard asymptotic behavior 

We briefly introduce the notion of standard asymptotic behavior (see 0). Speaking about a solution 
having standard asymptotic behavior, we mean first of all that this solution has the asymptotics 
and other properties that we present now. 

Fix Eq G C. Let D G C a, simply connected domain containing no branch points. Let /t be a branch of 
the complex momentum continuous in D and ip± the canonical Bloch solutions normalized in p^ G D. 
We say that a consistent solution / has standard behavior / ~ e^ '^^'^'>^'^ipj,.{x,p,E), respectively 

f ~ e~i E) in D if 
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• there exists a complex neighborhood Vo of Eq and X > 0 such that / is a consistent solution of 
equation (If.dll for any {x,ip,E) G [—X,X] x D x Vo; 

• for any x G [—X,X], the function {{(f,E) f{x,ip,E,e)) is analytic on D x Vo; 

• for any A, a sub-admissible domain of D, there is a neighborhood Va of E^ such that 

y{x,v.E)G[-X,X]xDxVA, f{x,^,E,e)=ei^^^^'^^^^{ilj+{x,v.E)+o{l)), e^O (5.10) 
respectively 

y{x,^,E)G[-X,X]xDxVA, /(a;,v?,£;,e) e^O (5.11) 

• the asymptotics (ICT71) and 15.1111 are uniform on [—X, X] x D x Va', 

• the asymptotics (ICTHi and imii can be differentiated once in x. 

5.5.2 

Let us formulate the Theorem WKB on a compact domain. 

Theorem 3. 

We assume that V satisfies (Hy) and that W satisfies (iLvv,r)- Eix X > 1 and Eq G C. Let 
K C Sy be a bounded canonical domain with respect to k. There exists eo > 0 and a consistent ba¬ 
sis {/+(x, X, e),/_ (cc, (/?, i?, e)} of solutions of (11.311 . having the standard behavior (15.101) et (15.1111 in 
K. 

For any fixed a: € M, the functions ip f±{x, ip, E, e) are analytic in K. 

5.6 The main geometric tools of the complex WKB method 

In this section, we introduce the main geometric tools of the complex WKB method. To do that, we 
recall some ideas of IH1IE1CIE3 

5.6.1 Stokes lines 

The definition of the Stokes lines is fairly standard, HE]. The integral ip i—> K{u)du has the same 
branch points as the complex momentum. Let ipQ be one of them. Consider the curves beginning at ipQ, 
and described by the equation 

Im f {k{^) - K{(po))d^ = 0 (5.12) 

Jipo 

These curves are the Stokes lines beginning at tpo. According to equation (ITnli and equation (14.1211 . the 
Stokes line definition is independent of the choice of the branch of k. 

Assume that W'{ipo) 0. Equation (I4.1)|l implies that there are exactly three Stokes lines beginning at 
ipQ. The angle between any two of them at this point is equal to 


5.7 Lines of Stokes type 

We recall that D C Sy is a simply connected domain containing no branch points. Let 7 C U be a 
smooth curve. We say that 7 is a line of Stokes type with respect to k if, along 7, we have 

/ rv \ / rv> \ 


either Im \^j K(u)duj = Const or Im {k{u) — TT)duj = Const 


5.8 Pre-canonical lines 

Let 7 C d be a vertical curve. We call 7 pre-canonical if it consists of union of bounded segments of 
canonical curves and/or lines of Stokes type. 
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5.9 Some branches of the complex momentum 

In this section, we describe different branches of n near the branch points described in l2.4l The geometrical 
configuration is similar to the one studied in [^. 

5.9.1 Different cases 

We assume that {Hw,r), {Hw,g), and {Hj) are satisfied. To study the geometrical tools of the WKB 
complex method, one needs to specify the properties of Im k and Re k. We know that K{(pf) = 0 [ 7 r] (see 
section^J. We consider two cases: either K{(pf) = 0[27r] or = 7r[27r]. 

We define S- the open domain delimited by the real line at the bottom and by E+ to the right: 

S- = {if — r ] if € r G R(i_} n Sy (5.13) 

Similarly, we define 5'+ the open domain delimited by the real line at the bottom and by E+ to the left: 

5'+ = {(/? + r ; If G r G R^} n Sy (5-14) 

The domains 5"+ and S- are shown in figure El 
We prove the following result. 

Lemma 5.2. There exists a branch m of the complex momentum such that 

1. Im Ki(ip) > 0 for ip G S-, Ki{ipf + *0) = 0 and Ki{ipi — 0) = tt, 
ou 

2. Im Ki{ip) < 0 for ip G S-, Ki{pf + iO) = tt and Ki{pi — 0) = 0. 

Proof 

• First, we specify the sign of Im Ki. The set {E — IF)(R — [pf, belongs to a gap G. We define 

A_ = {E — IF)(5'_). We prove that A_ is a connected domain which intersects with R only in the 
gap G. According to assumption {Hw,g) (12. 2. 311 . there exists a sequence of vertical curves Ej, such 
that: _ 

A_nR={E- W)((-oo,(^-] U [:^+, +oo)) UiE- W)(E+). 

{E — 1F)(E^) is a connected domain of R; it contains at least a point of G and does not intersect 
with da{Ho). Consequently, {E — IF)(E^) belongs to G and: 

A_ nR = G. 

We fix on A_ a continuous branch of the quasi momentum k. The sign of Im k does not change 
since {E — IF)(5'_) does not intersect with a{HQ). If we define Ki{p) = k{E — W(p)), Im we can 
assume that Im > 0 on S-. 

• Now, we consider Re Ki. According to sectional we can choose the branch Ki such that +i0) G 

{ 0 , 7 r}. First, we study the case Ki{pf) = 0; this assumption implies two possibilities. 

1. The point 0 is a minimum for W and the band B in (Hj) is in the form [£’ 4 ;+!, £ 4 ;_|_ 2 ]. The points 
Er and Ei satisfy Er = £ 4 ;+! and Ei = £ 4 z+ 2 . There exists a neighborhood V of [</?“,0]Ucr such 
that {E — W){S-nV) C C+\R. Actually, in the neighborhood of 0, we have Im (£ — IF(i^)) > 0; 
according to {Hw,g), there exists a neighborhood F of [(/?“, 0]Ucr such that [E — W){S-{^V) does 
not intersect R. By continuity of the mapping p 1 -^ Im (£ — W{p)), the sign of Im (£ — W{p)) 
remains positive on S- H V. There exists a branch k of the quasi-momentum such that 

Im k{£) > 0 for Im £ > 0 and fc(£„ -|- iO) = 0, k{Ep + iO) = tt. 

We define Ki{p) = k{E — W{p)). 

2. The point 0 is a maximum for W and the band B is in the form [£ 4 ;_|_ 3 , £ 4 ^+ 4 ]; the points Er and 
Ei satisfy Er = £ 4;+4 and Ei = £ 4 ^+ 3 . Let k be the branch of the quasi-momentum such that 
Im k{£) > 0 for Im £ < 0, and k{En) = 0; then k{Ep) = tt. We define Ki{p) = k{E — W{p)). 
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Figure 6: The domains S- and 5+ 

The case Ki((p~) = tt is similar. 

This completes the proof of the lemma. 

For the sake of clarity, for all the proofs, we shall consider the case: 

Ki{ip~ + iO) = 0 et Ki{ipi — 0) = TT (5.15) 

The arguments in the second case are similar and we will not give the details. 

5.9.2 Other branches of the complex momentum 

The properties of the complex momentum near the branch points tpi, Tpi, (pf- are determined by the 
behavior of k near Er and Ei. 

Now, we describe other branches of k, which are obtained from the branch Ki IT/eiTima, 15.211 by analytic 
continuation. We consider the case 

• We denote by Kg the continuation of Ki to the domain {Re (ip) < Re (p~)}. Kg satisfies 

Im {Kg{p)) > 0 for {Re (p) < Re 
Re {Kg){p) —*■ 0 as Re (i,^) —*■ —oo. 

Kg is the continuation of Ki through (—oo, p'^\. 

• We consider the strip Sy cut along (S\o-) U (S\ct) U (—oo, p~) U (i^+, +oo). We always denote by Ki 
the continuation of Ki through C. 

• On {Re {p) > Re ((/?+)}, we fix a continuous branch k^ with the conditions: 

Im {Kd{p)) > 0 for {Re {p) > Re (iy9+)} 

Re {Kd){p) —*■ 0 as Re (p) +oo. 

Kd is the continuation of Ki through 5+. 

Here, we describe the behavior of the different branches of k. 


Vp e 

Kg(p) = K,(p) ; 

VpeS-, 

1 

II 

c, 

-Ki{p). 

(5.16) 

Vp e S+, 

Kd(p) = -Ki(p) 


Kd{p) = 

Ki{p). 

(5.17) 


5.10 Stokes lines 

This section is devoted to the description of the Stokes lines under assumptions {Hw,g) and (Hj). We 
describe the Stokes lines beginning at p~, p^, pi and Since W is real on the real line, the set of the 
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Stokes lines is symmetric with respect to the real line. 

First, Ki is real on the interval C M; therefore, [ip ~is a part of a Stokes line beginning at 

ip~. The two other Stokes lines beginning at ip~ are symmetric with respect to the real line. We denote 
by b the Stokes line going upward and by b its symmetric. Similarly, we denote by a and a the two other 
Stokes lines beginning at ; a goes upwards. 

Consider the Stokes lines beginning at pi. The angles between the Stokes lines at this point are equal to 
27r/3. So, one of the Stokes lines is situated between E and e~E. It is locally going to the right of E; 
we denote by d this line. Similarly, we denote by e the Stokes line between E and e 5“ E. Finally, we 
denote by c the third Stokes line beginning at pp c is going upwards. 

By symmetry, we denote by c, d and e the Stokes lines beginning at 'pi. 

We describe the behavior of a, b, c, d and e in the strip Sy- We have represented these lines in figure [T] 
In this figure, we have precised the values of k in the branch points. 

Lemma 5.3. We assume that V, W and J satisfy (Hy), {Hy/) and {Hj). Then, the Stokes lines 
described in figure^have the following properties: 

1. a stays vertical; it intersects {Im (p) = Y}. 

2. b stays vertical; it intersects {Im {p) =Y}. 

3. d intersects a above pf; the segment between pi and this intersection with a is vertical. 

4-. e intersects b above pf; the segment between pi and this intersection with b is vertical. 

5. c stays vertical; it intersects {Im {p) = Y} and does not intersect a. 

6. a and c do not intersect one another in the strip Sy. 

7. b and c do not intersect one another in the strip Sy. 

Proof First, we note that a Stokes line can become horizontal only at a point where Im k = 0, 
i.e. at a point of the pre-image of a spectral band. Besides, a Stokes line beginning at pf (respectively 
at p 2 or is locally orthogonal to i k{p) (respectively i {tt — k{p))). 

We first prove 1). According {Hj), the pre-image of the spectrum is [pf ,pf]yja. So, a becomes horizontal 
only if it intersects a. Let us prove by contradiction that it is impossible. Let us assume that a intersects 
a in Pa, then: 

rVa. rVa, 

Im / K{u)du = 0 = Im / K,{u)du 

JJo, along cr 

(Re K{u))d{hn (u)) < —ki{E — tT_)Im pa < 0 

which is impossible. Therefore, a stays vertical. Moreover, as ^ oo, p € Sy, Im {iH) 0. Thus, a 
admits a vertical asymptote and intersects {Im {p) = Y}. 

Similarly, we prove 2). 

To prove 3), we consider the Stokes line d. If a and d do not intersect one another, then d intersect either 
a or [0, pf]. In this case, we denote by pd the intersection between d and cr and we have: 

rVi rVi 

Im / (k(m) — Tr)du = 0 = / Re (k(m) — 7r)d(Im (u)) < 0 

Jpd Jpd 

Consequently, d and a do not intersect one another. Before its intersection with a, d does not intersect 
the pre-image of a spectral band and it stays vertical. We prove similarly the properties of e. 

We prove now 5). c is going upwards, c does not intersect the pre-image of a spectral band in {Im p S 
]Im pi,Y\} and c stays vertical. 

We prove 6) by contradiction. Let us assume that there is (^a G a C c. Then, we compute: 

rVa r rVa 

Im / K{u)du = 0 = Im / K{u)du + Im / K{u)du 
Jo J <7 J ‘Pi 
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Figure 7: Stokes lines 


First, Im K{u)du = Trim {ipa — <pi) > 0 and Im n{u)du = Re K(u)d(Im u) > 0. 
which is impossible. So, a and c do not intersect one another in Sy- <) 

In the following, we choose Y g] sup Im ipi{E), y[. 

EeJ 


5.11 Construction of a consistent basis with standard behavior in the neigh¬ 
borhood of the cross 

In this section, we begin with constructing a canonical line near the cross. To do that, we follow the 
methods developed in |5]. 

5.11.1 General constructions 

We first recall some general geometric tools presented in |5], section 4.1. 

• We first introduce the idea of enclosing canonical domain. 

Definition 5.5. Let ^ CL D he a line canonical with respect to k. Denote its ends by and ^ 2 - Let 
a domain K G D be a canonieal domain corresponding to the triple k, and ^ 2 - If 1 K, then K 
is called a canonical domain enclosing 7. 

We have the following property: 

Lemma 5.4. m 

One can always construct a canonical domain enclosing any given compact canonical curve located 
in an arbitrarily small neighborhood of that curve. 

Such canonical domains, whose existence is established using this lemma are called local. 

• To construct a canonical domain, we need a canonical line to start with. To construct such a line, 
we first build pre-canonical lines made of some “elementary” curves. Let 7 C D be a vertical curve. 
We call 7 pre-canonical if it is a finite union of bounded segments of canonical lines and/or lines of 
Stokes type. The interest of pre-canonical curves is the following: 

Lemma 5.5. m 

Let j be a pre-canonical curve. Denote the ends of ^ by o-nd xit. Fix V G D, a neighborhood of 
7 and Va G D a neighborhood of ■ Then, there exists a canonical line 7 connecting the point to 
some point in V/. 

5.11.2 Constructing a canonical line near the cross 

Here, we mimic the construction of P], section 4.2. We assume that assumptions (i7y), {Hw,r), {IIw,g) 
and {H'j) are satisfied. We now explain the construction of a canonical line going from {Im ^ = —Y} to 
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Figure 8: A canonical curve 


{Im ^ = F}. First, we consider the curve f3 which is the union of the Stokes line b, the segment 0] 
of the real line, the closed curve cr+ and the Stokes line c. 

We now construct a a pre-canonical line close to the line f3. We prove: 

Proposition 2. Fix d > 0. In the 5-neighborhood of P, there exists a a pre-canonical line with respect 
to the branch k connecting to ^2 and having the following properties: 

• at its upper end, Im 1,2 = y, 

• at its lower end, Im = — Y, 

• it goes around the branch points of the complex momentum as the curve shown in figure\^ 

• it contains a canonical line which stays in S-, goes downward from a point in S- to the curve a and 
then continues along this curve until it intersects the real line. 

Proof The proof of this Proposition is completely similar to the proof of Proposition 4.2 in [^. It 
consists in breaking down a in “elementary” segments. We do not give the details. <0 
An immediate consequence of Proposition |21 is the following result: 

Proposition 3. In arbitrarily small neighborhood of the pre-canonical line a, there exists a canonical 
line 7 which has all the properties of the line a listed in Proposition\^ 


5.11.3 Some continuation tools 

In this section, we recall some continuation tools; these tools are developed in [ 3 . 

1. Now, we present the continuation lemma on compact domains. We recall that q is defined in 
Lemma 5.6. m Let if -, (/?+, ipQ be fixed points such that 

• Im if- = Im (f^; 

• there is no branch point of ip n{ip) on the interval 

• (^0 e {F-F+),q{Fo) 0. 

Fix a continuous branch of k on [ip-,ip+\. Let f{x,ip,E,e), f±{x,p, E,e) be solutions of (11.311 for 
if S [if -, ipjf\ and x € {—X, X\ satisfying (11.511 and such that: 

(a) f{x,ip, E, e) = e^-^'co p,E) + o(l)) pour p G [p-,pq] for p G [p-,pq] when e ^ 0 

and the asymptotic is differentiable in x; 

(b) f±{x,p,E,e) = o{\)) for p G [p-,pff\ when e ^ 0, and the asymp¬ 

totic is differentiable in x. 

Here, il)± are canonical Bloch solutions associated to the complex momentum k. 

Then, 
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Figure 9: Stokes lines and Stokes Lemma 


• if Im {k{ip)) 


• if Im (k(<p)) 


> 0 for all ip € [p-, P+\, there exists C > 0 such that, for e > 0 small enough, 


dx 


{x,p,E,£) 


+ \fix,p,E,e)\ < Ce^ 1 -^"* >^<,n)\du^ ^ ^ [po,p+\] 


< 0 for all p G [p-, P+\, then 

f{x,p,E,£) = +o(l)), p G [pq,P+\, 


(5.18) 

(5.19) 


and the asymptotic is differentiable in x. 

Intuitively, this lemma means that a function / has the standard behavior along a horizontal line as 
long as the leading term of its asymptotics is growing along that line. For analogous results with real 
WKB method, we refer to m- 

2. The estimate we obtained in Lemma 15.61 can be far from optimal. The Adjacent Canonical Domain 
Principle gives a more precise result: 

Proposition 4. e Assume that a solution f has standard behavior in either the left hand side or 
the right hand side of a constant neighborhood of a vertical curve 7. Assume that 7 is canonical with 
respect to some branch of the complex momentum. Then f has standard behavior in any bounded 
canonical domain enclosing 7. 

3. The last tool we shall need in the sequel is the Stokes Lemma. 

Notations and assumptions: 

Assume that ^0 is a branch point of the complex momentum such that kF^(Co) 7 ^ 0. There are three 
Stokes lines beginning at ^o- The angles between them at ^0 are equal to 27r/3. We denote these 
lines by cri, (J2 and 0 - 3 , so that cri is vertical at ^o- Let P be a neighborhood of ^ 0 ; assume that V is 
so small that cri, (72 and (73 divide it into three sectors. We denote them by Si, S2 and S3 so that 
Si be situated between cri and (72, and the sector S 2 be between (72 and (73 (see figure 0 . 

We recall now the result: 

Lemma 5.7. Let V be sufficiently small. Let f be a solution that has standard behavior f = 
e 7 1“^ '^("i‘^“('!/;_|_(a;, p) + o(l)) inside the sector Si U (72 U S 2 of V. Moreover, assume that, in Si near 
(7i, one has Im k > 0 if Si is to the left of ai and Im k < 0 otherwise. Then, f has standard 
behavior f = ”*^“i'^“(^+(a;, (/?) + o(l)) inside V\ai, the asymptotics being obtained by analytic 

continuation from Si U (72 U S2 to V\ai. 


5.11.4 Construction of a basis with standard asymptotic behavior near the cross 

We prove the existence of a consistent basis with standard asymptotic behavior near the canonical line 
a. Let a be the curve described in Proposition |31 According to Lemma 101 we can construct a local 
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canonical domain Ki enclosing a. 

Proposition 5. Assume that {Hy), {Hy/,r), {Hw,g) cLi^d, (Hj) are satisfied. Fix Eq € J, X > I and 
Y €]0, Y[. Then, there exist a complex neighborhood Uq of Eq, a real number eo > 0 and a function fi 
satisfying the following properties: 

• The function {x,(p,E,e) ^ fi{x,tp,E,£) is defined on R x Sy x £o[- 

• For any a; e K, for any £ g] 0, eo[) the function {{(p, E) i—> fi(x, tp, E, e)) is analytic on Sy x Uq. 

• For (x,p,E) G [—X, X] X Ki X Ido, the function fi has the asymptotic behavior: 

fi{x,p,E,£)=eihf<'^)'^-^{pj}^[x,p,E)+o{l)), £^0. (5.20) 

• The asymptotics (lonii are uniform in {x,p,E) € [—X,X] x Ki x Uq. 

• The asymptotics can be differentiated once in x. 

• There exists a real number Ui S { — 1,1} such that the function fi satisfies the relation: 

wifijf) = w{fi{-,p,E,£),fi{-,lp,E,£)) = ai{k[wi){E - W {Q)) 

The end of this section is devoted to the proof of Proposition [3 This Proposition mainly follows 
from Theorem 121 

5.11.5 Existence of fi 

The domain Ki is a local canonical domain. According to Theorem |21 we can build a function fi such 
that, on Ki, fi has the following asymptotic behavior: 


Let us normalize /i in 0. 

5.11.6 Computation of the Wronskian w{fi,f*) 

To finish the proof of Proposition |21 it remains to compute w{fi, /*). 

Let R+ be a small enough rectangle to the left of a+, so that R+ C Ki H Sy. We define R = R+ U i?_; 
we study the behavior of fi and (fi)* in R. 

• First, by construction, in Ry, the function fi satisfies: 

• To the right of a-, the function fi satisfies 

with Im Ki < 0. According to Lemma l5.6l we know that, in S-, the function fi admits the asymptotic 
behavior: 

• Thus, the function fi has the standard asymptotic behavior in R. 

• Now, we study the behavior of /*. To do that, we start with describing the main objects related to 
Ki in R. Let ki be the branch of the quasi-moment of 63, analytically continued through [Er,Ei\ 
and satisfying: 

ki{Er) = 0 and ki{Ei) = tt 
ki is real on [Er,Ei]. Therefore, ki satisfies: 

h(£) = k^. 
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The branch Ki satisfies Ki{(p) = ki{E — W{(p)). The associated canonical Bloch solutions are such 
that: _ 

Therefore, we have in R: 

e R. (5.21) 

Besides, since k[ is real on the band, there exists a real number Ui S {—1,1} such that: 

(5-22) 

We shall precise this coefficient in section En 
We compute: 

w{r+{-, E, e), E, e)) = qi{0)q*0),'^^_{-,0))g{(p, E, e). 

Since 0), '!'!_(•, 0)) = —rc('I'!|_(-, 0), 0)), the term g{(f, E, e) satisfies: 

g*{(p, E, e) = g{(f, E, e) = g{ip, E, e), 
g{ip,E,e) = [1 + o(l)]. 

Since the Wronskian is analytic and £-periodic, this asymptotic is valid in Sy- 

Since g* = g and g = [1 + o(l)], there exists an analytic function {(p, E) h{(p, E, e) on Sy such 

that: 

- g{^,E,e) = h{(f,E,e)h*{(p,E,e), 

- H‘P,E,e) = [l + o(l)]. 

We slightly deform fi, i.e., we replace fi by ; the basis {fi, f*} is consistent. 

This ends the proof of Proposition E 


6 Consistent Jost solutions of (II.HI) 

This section is devoted to the proof of the following result. 

Proposition 6. We assume that {Hy), {Hw,r) and (Hj) are satisfied. Fix X > 1 and A > 1. Then, 
there exist a complex neighborhood V = V of J, a real Sq > 0, a constant C > 0, two complex numbers 
mg, md and two functions {x,(p,E,e) hf{x,ip,E,£), {x,ip,E,£) hy{x,(p,E,£) such that, if we 

define 

S Sy Re ifi < —et € Sy ; Re p > C£~^^'^ , 

then 

• The functions {x,p,E,£) i—> hf{x,p,E,£) and {x,p,E,£) i—> hy{x,p,E,£) are clearly defined and 
consistent on M x Sy x Vx]0,eo[- 

• For any X G [—X, X] and e g] 0, £o[, {p,E) hf{x,p,E,£) and {p, E) hy{x,p, E,£) are analytic 

on Sy X V. 

• The function X hf(x,p,E,£) (resp. x i—> hy{x,p, E,£)) is a basis of J- (resp. J+). 

• The functions hf and hy have the following asymptotic behavior: 

hf{x, p, E, £) = e~ ■^”*9 p, E){1 + Rg{x, p, E, £)), (6.1) 

and 

hyix, p, E, £) = e" '^^"^‘'“'(/'+(a:, p, E){1 + Rd{x, p, E, £)), (6.2) 

where 
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Rg and Rd satisfy: 


3M > 0, 
3M > 0, 


Ver g]0, £o[, 

V£ G]0, £o[j 


,Va;e[-X,X], yEeV, y>peBI, 
,Vxe[-X,X], yEeV, 


I / M ^ 

M 


— The functions '0+ and tp- are the Bloch canonical solutions of the periodic equation eu defined 
in section o 

• The asymptotics E2P and 16. iil) may be differentiated once in x. 

• There exist two real numbers Ug G { — 1,1}, cfd G { — 1,1}, an integer p and two functions E otg[E) 
and E i—*■ ad{E) such that: 

1. For any e g]0,£o[, a; G K, i? G V,et ip G Sy ,we have: 


ag{E)hf {x,p,E,e) 


= i<Tge ^‘^^'^^ag{E)hf{x,p,E,e) 


(6.3) 


ad{E)h'l{x,p,E,£) = iadei^^'^'^adiE)h'l{x,p, E,e) (6.4) 

2. The functions ag and ad are analytic and given by (lO^ and (i?r^ . They do not vanish on V. 

We shall construct some consistent Jost solutions of (tot . To do that, we regard equation as 
a perturbation of equation mi with £ = E. We adapt the construction of Jost functions developed in 
nniini. Precisely, we look for solutions of ED in the form : 

p9 =e-*'=(^)‘^/"0O(a:,£;)(l + o(l)), x ^-oo, 
pd ^ (x, E){1 + o(l)), a; ^ +oo. 

Since the functions {x,p,E,e) i?) are consistent, they allow us to construct a consis¬ 

tent resolvent for the periodic equation. Using this property and the fact that equation is invariant 
by the consistency transformation (a:, p) i—>- {x— 1, p + e), we obtain the consistency of the Jost functions. 


6.1 Construction of the Jost functions 

We start with constructing Ef. The construction of Ef is similar. Since the parameter E lies in the 
neighborhood of a gap, Im k{E) is non zero; the function Ff is therefore exponentially decreasing and 
goes to zero as x goes to — oo. Such a solution is called recessive. 


6 . 1.1 

On a small enough complex neighborhood of J, V = V, one can fix a determination k of the quasi¬ 
momentum such that: 

Im k{E) >P>0, VU G V. 

Fix mg in Sy such that: 

• The point mg is not a branch point of k. 

• It satisfies Im mg > 0, k'^{mg) 0. 

• The domain {p G Sy ; Re ((^ — mg) < 0 and Im (p — mg) > 0} does not contain any branch point 
of K. 
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We define Eg = E — W{mg). We denote by tjj^ the analytic Bloch solutions of equation 12.111 normalized 
at the point Eg ( k'{Eg)) ^ 0). These solutions are constructed in Lemma [4.11 

with pI{x+ 1 ,E) =pI{x,E). 


We define : _ 

tlj±{x,ip,E,e) = e^^’^^^'>^^+^'>pl{x,E) = 

We consider the resolvent R oi Hq ; 

'tpl{x,E)'tlj°_{x',E)-'ipl{x',E)tlj°_{x,E) , 


{Rg){x) = - f 
J — ( 


{k'wo){Eg) 


-g{x')dx 


6.2 

Since '0- goes to zero as x goes to — oo, we look for a recessive consistent solution / of 1E3 in the form 

/(x, ip, E, e) = 'll:- [x, ip, E, e) + R[W{ex + ip)fix, ip, E, e)]. (6.5) 

We define f{x,ip,E,e) = if, e); equation 16.511 is transformed into: 


f{x,ip,E,e)=p^_{x,E)+ f A{x,x',E)W{ex'+ ip)f{x',ip,E,£)dx' 

*/ —OO 


where the function A satisfies: 


A{x, x' ,E) = 


^ 2 ik{E){x X )pO_,^{x,E)p'^_{x' ,E) — p'i^{x', E)p'^_{x, E) 


ik'wo){Eg) 

Since Im k{E) > /3 > 0 for E G V, there exists a constant C > 0 such that: 

Vx > x', WEgV, |A(x, x', if)| < C 


( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


6.2.1 

Fix Xq G M and a > 0. If / is a real interval, we define: 

Rf = {ip G Sv ; Re ip G /}. 

Let i?((— oOjXq] X R[-a,a]) the set of bounded functions {/ : (x,:^) f(x,ip)} on (—oo,Xo] x R[-a,a]- 
The set ii((—oo,Xo] x R[_a^a]) equipped with the norm 

||/||oo= sup lf(x,‘p)l 

xC (—cxd.Xq] ,R.e tpe[—a,a] 

is a Banach space. 

We define the integral operator Tb by: 

Tb : B((-00, Xo] X R[_a^a]) —> -B((-00,Xo] X i?[_o_o]) 

f ^ F 


where F{x,p) = f A{x,x ,E)W{£x + ip)f{x ,ip)dx . (6.9) 

J —OO 

The operator Tb is a bounded operator on ii((—oo,Xo] x R[_a^a]) and satisfies the estimate: 

Vx G (-oo,Xo], V(^ G i?[_a,a], (/) (a^^, V^)I < C'll/II oo f | IT(£:x' + Re (i^) + zim ((^))| dx'. 

J —OO 


\\TE{f)\\oo < 


M 


snp I / M 

e ,e(_oo.Xo],Re ^e[-a,a] + Re (<^)| 


S-1 ' 
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6.2.2 


Fix A > 1. There exists a constant C > 0 such that: 

|Xo| >Ce-^ ^ lllTslll 

We rewrite 116.611 in the form: 

il-TE)f=p°_ix,E) 

The operator is then invertible on B{{—oo,Xo] x We define: 

Fl(x,(p,E,e) = (1 - TE)~^p°_ix,ip,E,e) 

We now give some properties of F"®. 


6.3 Properties of F® 

6.3.1 Asymptotic behavior in x 

Substituting 

(1 - Te)-^ = 1 - (1 - Te)-^Te 

in equation (EUl), we obtain: 

El {x, p, E, e) = £;)(! + Rg{x, p, E, e)), 

with 

M 

\Rg{x,p,E,e)\ < 

for x G (—oo, Aq] and p € R[-a,a]- 

The function F® is therefore in the dost subspace JL of equation dOJ. 


( 6 . 10 ) 

( 6 . 11 ) 


( 6 . 12 ) 


6.3.2 Study of the consistency 

We assume that a > 1 and e < 1. We now prove that the function F® is consistent. 

We denote by G the function: 

G : (x, (p, E, e) 1 -^ G(x, (p, E, e) = Ef_(x + 1, (p — e, E, e) 

G is defined for x G (—oo, Aq —1] and ip G F[_a+i_a_i]. Moreover, the function G belongs to F((—oo, Aq — 
1] X F[_(j+i_(j-i])- 
We define the operator: 

Te '■ F((-00, Aq - 1] X F[_a+i^a_l]) ^ F((-00, Aq - 1] X F[_a+l,a- 1] ) 

f ^ F 


where F(x, (p) = f A{x,x', E)W{ex'+ (p)f{x', ip)dx'. (6.13) 

J —OO 

Since F((—oo,Ao] x F[_o_o]) C F((—oo,Ao — 1] x F[_a+i_o_i]) and according to equations 16.911 and 

(ICTill . the operator Te is an extension of the operator Te- Let us denote by F® the restriction of F® to 

( 00, Aq 1] X 

We compute in F((—oo, Aq — 1] x F[_a+i,a-i]): 

{Te{G)){x, ip, F, e) = (rE(F® ))(a: + 1, (/?-£, F, e) 
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This leads to: 

{{l-TE){G)){x + l,ip-e,E,e) = {{I- TE){Fi)){x + - e,E,e) =p°_{x + l,E) =p°_{x,E), 

The functions and G satisfy the relation: 

((i-r^)(G)) = ((i-i^)(T^)). 

For a sufhciently small eo, the operator Te satisfies, for any e g]0,£o[: 

The operator (1 — Te) is invertible in i?((—oo, Xq — 1] x and: 

F^=G. 

For If G R[-a+i,a-i]^ the functions F® and G coincide on (— 00 , Xq — 1]; according to the Cauchy-Lipschitz 
Theorem, they coincide for a: G K. Fix x G M; F® and G coincide for ip G F[_a+i,a-i]- By analyticity, 
they are equal for ip G Sy- 


6.3.3 Asymptotic behavior in ip 

We use now the consistency of F® to compute its asymptotics as Re ip goes to — 00 . Fix X > 0. We 
study F® for x G [—X, X]. The function F® is consistent, and: 

F®(x,(p,F,e) = F®(x+ [Re {ip)],E,e) 






ejex + [Re ' 

Re < -Ce-^ ^ F® (x, p, E, e) = e){1 + F^(x, ip, E, e)), 

\Rg{x,ip,E,e)\ < 


As a result, there exists a constant C such that: 

X 


where 


M 
e|Re ip\ 


s-i ■ 


for X G [—X, X] and Re (/? < —Ge . 

We define = {ip ^ Sy ■, Re 

6.4 Renormalization of Ff 

We now renormalize F®. We define: 

f^{x,ip,E,e) = e~~^ (6.14) 

where the integral °° [«^(u) — k{E)\du is taken in the upper half plane. The function F 1 -^ J^°°[K — k{E)] 

is analytic on V. For </? G Ff, we have: 

f^{x,ip,E,e) = e~‘ + o{l)) 


Since the function is analytic and since W{ip) = 0(£=-i) for ip G F|, we get: 

yp G Ff, i^-{x, P,E) = ip°_{x,E - W{p)) = tjj°_{x,E){l + 0 ( 1 )) 
We finally obtain that, for x G [—X, X] and p G B§: 

f-{x, p, E, e) =e~~ F (x^p^E)(l + o(l)) 


(6.15) 
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6.4.1 Symmetries 

Let 7 be a complex path and / be an analytic function on 7 . We have: 


f fiz)dz= f_f*{z)dz. 
J -y J ^ 


(6.16) 


Since J satisfies according to equation 14.1211 . there exists an integer p such that: 


k{E) + k*[E)=2pT:. (6.17) 

We recall that the functions iu± associated to k are defined by equation We consider a path 7 ^ 

such that: 

• The path 7 ^ connects rUg to mg and is symmetric with respect to the real axis. 

• The path 7 ^ does not contain any branch point of k and any pole of U!±. 

We fix a continuous determination qg of ^/k^ on 7 g. According to relation 14.1211 . we have {k*)' = —k', 
which implies that there exists (Jg G { — 1 , 1 } such that: 


dg = 

The functions 'tp±(x, p, E, mg) satisfy the relation: 

'tp±(x,(p,E,mg) = itTge^^^e-^^o p, E, mg). 

Besides, equations (IHT^ and BTtIi lead to the following relations: 

= — UJ+ ; / u!- = — I UJ-. 


la 


'la 


(6.18) 

(6.19) 

( 6 . 20 ) 


la 


'la 


According to {Hw,r), W* = W. By using 16.711 . we compute: 


A{x, x', E) = A{x, x', E). 

The operator Te satisfies: 

TEin = [TEif)]*. 

Consequently, according to (OUll and (Einii, we obtain that, for if in V, a: in R and (p in B®, 
(F»)*(x, p, E, e) = F^{x,p,E, e) = (x, E, e). 

This leads to: 


( 6 . 21 ) 


where 


AE) 


^ " a*g{E) - 

^ p-4i{S:^giAA-P'^)<lu+p'K{mg+m^)) . 


( 6 . 22 ) 


Similarly, we fix m^ in Sy such that: 

• The point md is not a branch point of k. 

• It satisfies Im md > 0, k'^{md) ^ 0. 

• The domain {ip G Sy ; Re ((/? — md) > 0 and Im {p — md) > 0} does not contain any branch point 
of K. 

We define Ed = E — W{md) and we define the function by: 

h%{x,p,E,e) F^(x,(p,A,e) (6.23) 

where the integral /^J°[«:(u) — k{E)\du is taken in the upper half plane. 

We consider the path 7 ^ such that: 
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• The path 7 ^ connects rrid to rrid and is symmetric with respect to the real axis. 

• The path 7 ^ does not contain any branch point of k and any pole of u}±. 

We fix a continuous branch qj, of on 7 ^. There exists a real number ad such that: 

Qd = ( 6 . 24 ) 


The function h'^ satisfies: 


where 


{hD* = iadei^P^^ 


oid{E) 


h 


d 

+ ’ 




We define the transmission coefficient: 


( 6 . 25 ) 


d{Lp,E,e) = w{aghf{-,ip,E,e),adh‘l{-,ip,E,e)) ( 6 . 26 ) 

We immediately deduce from Proposition ^ and Propositional that the eigenvalues of Hip g are charac¬ 
terized by: 

d{ip,E,e)=0 ( 6 . 27 ) 


6.5 Some remarks 

6.5.1 

The assumption {E[w,r) is not optimal. Actually, it suffices to assume that W is analytic real in Sy and 
that there exists a function / e L^(IR) such that : 

Vx € M sup \W{x + iy)\ < f (x). 
vel-Y,Y] 


6.5.2 

In equations JOl and llOll . we could have included the numbers iag and iad into the functions ag and 
ad, but we prefer showing the relations between qg and q*, qd and q^. 

6.5.3 

Note that this construction differs from the constructions of canonical domains in |^. Indeed, the domains 
on which we construct these functions depend on e. We shall extend these asymptotics on a fixed strip 
in the neighborhood of the real line (section EJ. 

7 WKB Theorem on non compact domains 

In this section, we prove a continuation result on non compact domains of Sy. This result is a general¬ 
ization on non compact domains of the method developed in and particularly of lyemma, 15.61 
We prove that the continuation of asymptotics stay valid on some half-strips {v? € Sy ; |Re > A}. 
To do that, we cover these domains by a countable union of small local overlapping canonical domains, 
called 5-chain (see section f?.1.511 . 

This principle follows the recent developments and improvements of the WKB method (see [7]). The 
idea is to get over the local notion of canonical domain in favor of maximal domains. These domains, 
constructed as union of local canonical domains are some domains on which a function keeps the standard 
behavior (see [Jj). 
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7.1 Continuation Theorem on non compact domains 

7.1.1 The main result 

We shall prove the following result: 

Theorem 4. Continuation Theorem on non eompact domains. 

Fix Y €]0, T[. Assume that V satisfies {Hy), that W satisfies {H\y,r) md that J satisfies (Hj). Then, 
there exist a real eo > 0, a eomplex neighborhood V of J and two real numbers Ag and Ad such that, if f 
has the following properties: 

• The function f(-,ip,E,e) is a consistent solution of (II. IH . 

• The function {p,E) i—> f(x,(p,E,e) is analytic on Sy x V for any x G [—X,X] and any e G]0,eo[. 
Then, 

1. There exists k a continuous branch on {Re ip < Ag} such that Im k > 0. Moreover, for any C < 
B < Ag, if the function f satisfies the asymptotic behavior 

f{x, if, E, e)=e~i tp, E) + rc{x, ip, E, e)) (7.1) 

with lira sup Taduy:{\rc{x,ip,E,e)\,\dxrc{x,ip,E,e)\} = Q, 

then, this behavior stays valid until B. Precisely: 

f{x,ip,E,e) = e" (V'- (x, £;) +rB(a;,v?,£^,e)) (7.2) 

with lim sup Tnax.{\rB{x,ip,E,£)\,\dxrB{x,ip,E,e)\} = Q. 

[-X,X]xfl(_o„,B]XV 

2. There exists k a continuous branch on {Re ip > Ad} such that Im k > 0. Moreover, for any C > 
B > Ad, if f satisfies the asymptotic behavior 

f{x,ip,E,e) = ei-^'"'^^'^^^'^{ip+{x,ip,E) + rc(x, p, E,£)) (7.3) 

with lim sup max{lrc(x, p, E, e)!, jd^rcix, p, E, e)!} = 0, 

<^^0 [-X,X]xfl[c,+oo)XV 

then this behavior stays valid until B. Precisely: 

fix, p, E, £) = eH" (^+(a:, E) + rB(x, p, E, e)) (7.4) 

with lim sup max.{\rBix,p,E,£)\,\dxrBix,p,E,£)\} = 0. 

<^^0 [-X,X]xfl[B,+oo)XV 

Theorem 0] and Proposition clearly imply Theorem [3 

7.1.2 Some remarks 

We shall prove Theorem 0] as W satisfies the weaker assumptions: 

(HI) W is an analytic real function in Sy- 
(H2)3C>0, 3 s > 1 such that V z G Sy, |W'(z)| < 

(H3) 3 f G L^(R) such that Vx G R sup |W(x + iu)| < f(x) 

ye[-Y,Y] 

The following lemma relates {Hw,r) and (771), (772) and (773): 

Lemma 7.1. Let W satisfy (Hwr) on Sy. FixY g]0, F). Then W satisfies (HI), (772) and (773) on 

Sy. 

Proof Assume that W satisfy (Hw,r) on Sy. We prove that W satisfies (772) on Sy by using the 
following lemma: 
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C > 0. 


Lemma 7.2. Let f be an analytic function on Sy such that |/(z)| < 

Fix 77 > 0. Then, 

VpeN* 3Cp>0/ yzGSy-r, 

1 + 2 :' 


Proof 

This result is a consequence of the Cauchy formula. We do not give the details. <C> 

- Clearly, W satisfies {HI) on Sy. 

- W satisfies {H3) with f {x) = ■ 

This completes the proof of Lemma 17. II <J> 


7.1.3 

Let us briefly outline the ideas of the proof. We shall concentrate on Bg = {ip G Sy; Re (ip) < Ag}. 
There are three steps. 

First we cover Bg with an union of overlapping local compact canonical domains 

In each canonical domain RTmj we can construct a consistent local basis thanks to Theorem El To 
compute the connection between the consistent bases of and Km+n, it suffices to do the product of 
the n transfer matrices between the canonical bases of two successive domains. The accuracy of the rest 
cannot be better than the sum of the accuracies obtained on each domain. Theorem El gives an estimate 
in 0 ( 1 ); this accuracy is insufficient when n goes to infinity. 

A refinement of the calculation of asymptotics in Theorem El is therefore necessary. We prove it by using 
the integrability of W. 

7.1.4 Branch points 

The following result specifies the location of the branch points of k. We recall that T{E) is defined in 

Lemma 7.3. Let V be a complex neighborhood of the interval J. Assume that W satisfies 

lim sup \W{x + iy)\ = Q, 

X^ + 00 


then: 

3A > 0 such that V A G V, p G T{E) n Sy |i?e ((/j)| < A. 

Proof Since V n da{Ho) = 0, there exists a > 0 such that: 

yEeV, ypGW, \E-Ep\>a. 

If Pp{E) satisfies E — W{pp{E)) = Ep, we get: 

yEGV, ypGN*, \W{pp{E))\>a. 

Finally, {u G Sy ; |IF(u)| > a} is a subset of a compact of Sy. This completes the proof of Lemma o 

❖ 

7.1.5 Uniform asymptotics on a d-chain 

First, we introduce a new definition. We remind that the width of a complex subset is defined in 15.911 . 
Definition 7.1. 5-chain of strictly canonical domains 

Fix Y g]0,U[. Fix E. Let D be a simply connected domain of Sy containing no branch points of the 
complex momentum. We fix on D a continuous branch k of the complex momentum. Let be a 

sequence of real numbers and K be a compact of Sy- 

[K + is called a 5-chain for E, n and D if it satisfies the following properties: 
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CXD 

1. U {K + Tr,)=D. 

n—0 

2. 3t > 0 such that Vn G N 1{{K + r„) n (i^ + Tn+i), Y) > r. 

3. The domain K is an union of curves 7 such that, for any n, 'y + Tn is a S-strictly canonical curve 
for K. 

K is called the fundamental domain of the 5-chain. Now, we have the intermediate result: 

Proposition 7. Assume that V satisfies (Hy) and that W satisfies {Hi), {H 2 ) and {H 3 ). FixY g]0, y[. 
Let V a complex neighborhood of J and D C Sy a domain with the following properties: 

• inf dist{D,ipr){E)} > C, 
peN*,Eev 

• there exists such that, for any E gV, {K -|- is a 5-chain for E and D. 

Fix (fio G D. 

Then, there exists £0 > 0 such that, for any n G N, there exist two functions {x,ip, E,e) 1 -^ 'if’^{x,g},E,e) 
with the following properties: 

• The functions {x,(p,E,e) ip{f{x, (p, E, e) are defined on R x {K -|- r„) x Vx]0,£o[ and form a 
consistent basis. 

• for any fixed a; G R, e G]0, £o[, the functions {(p, E) 1 -^ p, E, e) are analytic on {K -\- t„) x V. 

• for X G [—X,X], p G {K -|- t„) and E gV, the functions have the asymptotic behavior: 


'il:l{x,p,E,e) (^fj±{x,p, E) -G ^ ’ 

• The asymptotics are uniform in x, t, p G K -G Tn et E G V. 

• The asymptotics can be differentiated once in x. 


(7.5) 


The proof of Proposition [7| mimics this of Theorem 1.1 in |^. We omit the details and we refer to 
[H], section 4 for an analogous statement. 


7.2 Construction of a (5-chain of strictly canonical domains 

In this section, we shall construct a 5 -chain under assumptions {Hi), {H2) and {H3). 

Proposition 8 . FixY g]0,T[. Assume thatV satisfies {Hy), thatW satisfies {Hi), {H2) and {H3) and 
that J satisfies {Hj). Then, there exist a complex neighborhood V of J, two real numbers {Ag,Ad) G R^, 
a domain K C Sy and two real sequences such that: 

• for any E G V, there exists a continuous branch k on {p G Sy ; Re p G (— 00, Ag]} (resp. on 
{p G Sy ■, Re p G [Ad, oo)}j, 

• for any E G V, {K -G (resp. {K -G r^InGN/^ is a 5-chain for n, E and {Re p G (—00, 

(resp. [Re p G [^d,+oo)^/ 

The rest of the section insi is devoted to the proof of Proposition |S 1 This proof is based on 
elementary geometrical arguments. We prove the construction for Re G (—oo,^^]. 

7.2.1 Construction of 5-strictly canonical straight-lines 

We have defined the canonical lines in section 15 - 1 -dl and described them in terms of the vector t{p). 

We set a = i inf Im k{E) and m = 2 sup |Re k{E)\. 

E&J EeJ 

Since the mapping {E,p) ^ E — W{p) is continuous and since W{p) goes to zero when Re p goes to 
infinity, there exist a complex neighborhood V of J and a real number Ag such that: 

Vi? G V, V (/3 G (—00, Ag], Ke k{E— W{p)) G[—m,m], Im. k{E — W{p)) > a 

We set Bg = (—00, Ag] -I- i[—Y,Y]. The canonical curves for Re p in the neighborhood of —00 are 
described by: 
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Figure 10: The fundamental domain K 


Lemma 7.4. There exists Oq s]0, 7 r/ 2 [ such that, if "f is a smooth curve in Bg satisfying: 


G 7, a.Tg[t{ip)] g]6»o, 7 r /2 - 6»o[. 


(7.6) 


then, "f is a canonical line for k. 

Proof For arg(M) = 9 and cot0 g] — Tr+m-s have: 

Im ((k — S)u) >0 et Im ((tt — k + S)u) > 0 
Consequently, cot 9o = implies that 17.611 is satisfied. <0 

7.2.2 The fundamental domain K 

Let .fi = —iY and ^2 = iY. We denote by K the lozenge bounded by the straight lines containing and 
^2 whose guiding vectors have the affixes and 

We set [—Mo,uo] = K r\ {y = 0}. K is shown in figure [TT)I Fix x such that K + x C Bg] we shall show 
that iG + a: is a (5-strictly canonical domain. According to Lemma 17.41 it suffices to write K as an union 
of smooth curves satisfying (TTini . 

For any u G K, we consider a vertical segment containing u and included in K (see figure EJ. The 
broken line U U [^,^ 2 ] satisfies 17.611 . The relation is stable under small C^-perturbation; 

we slightly deform the line [^ 1 ,^] U U [^,.^ 2 ] to get a smooth curve which satisfies 17.611 . 
Consequently, K satisfies the following properties: 

- K n Sy contains a rectangle of width dry > 0. 

- 1{{K — nr]) n {K — (n + l)ry), Y) > ry. 

- K is the union of curves 7 such that 7 — nry is (5-strictly canonical for any sufficiently large n. 

7.2.3 Conclusion 

To finish the proof, it suffices to adapt the proof of Lemmain section 5.9 of |S], by using Proposition 
Inland Proposition [7| The convergence of the series of general term |3 replaces the compactness. We 
do not give the details. 
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8 Transmission coefficient. Equation for eigenvalues 

In Theorem 121 we have constructed two functions h^_ and h'^. We have defined the transmission coeffi¬ 
cient d{E, ip, e). We choose nig = —0 -I- iO and = 0 -I- iO. 

In Proposition 121 we have introduced a consistent basis (/i,/*) near the cross. To compute d{E,(p,e), 
we shall project the functions h^_ and h'i_ onto the basis (/i,/*)- 


8.1 Preliminaries 

8.1.1 Introduction. Notations 

Fix Y < Y and Eq G J . We have described in section [51 the complex momentum k and the related 
geometric objects. We recall that we consider the case We use the notations introduced in section 

121 The branch points are called ipf and ipi, T^. We have described the Stokes lines in section 15.101 
We have described in sections Ifi.l). 21 and IS. 1). II the different branches Ki, Kg and Kd- The branch Kg, resp. 
Kd, is defined and continuous on the domain {if G Sy', Re (/J < (^“}, resp. {Re </? > ip ^}. The branch Ki is 
defined and continuous on a neighborhood of the cross. The domain {E—W) {{ip G Sy, Re < ip~}) is a 
simply connected domain which intersects with real axis in only one gap. Thus, we can fix a determination 
kg of the quasi-momentum such that: 


kg{E-W{ip)) = Kg{ip). 

Similarly, we fix the branches ki and kd of the quasi-momentum such that : 

ki{E -W{ip)) = Ki{ip), kd{E -W{ip)) = Kd{ip). 


Finally, we set: 

q,{ip) = ^k[{E-W{ip)), qg{ip) = k'g{E - W{ip)) qd{ip) = k'^{E - W{ip)). 

Let G R such that ipg < ip~ and such that the interval [ipg,ip~] does not contain any pole of We 
define the path 7 g in the complex plane by: 

7g = [—0 -I- iO, ipg + iO] U [ipg — iO, —0 — iO]. 

Similarly, fix G R such that ipd > ip'^ and such that the interval [ip{:, pd\ does not contain any pole of 
uj±. We define the path ^d in the complex plane by: 

7<i = [0 + iO, ipd + iO] U [ipd - iO, 0 - *0]. 

In the following section, we explain the choice of the determinations qi, qg and qd- 

8.1.2 The determination q 

We recall that there exists a real number ai G { — 1,1} such that: 

— (8.1) 
di 

We refer to section mm 

The Wronskian satisfies w{fi, {fi)*) = ai{wok{){E — kF(0)). 

The number ai depends on the sign of k' along the band B: 

• If the band B can be written [i? 4 p+i, T^ 4 p+ 2 ], then k' > 0 on B and ai = 1. 

• If the band B can be written [i54p4-3, i? 4 p+ 4 ], then k' < 0 on B and di = — 1. 
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We fix the branch qg such that qg = qi in S- and such that qg is analytically continued in S Sy', Re q} < 
According to relation the branch qg satisfies: 

(fg = mqg ( 8 . 2 ) 

Similarly, we fix qd such that qd = qi in Sy and such that qd is analytically continued in {(/? S Sy; Re > 
}. The branch qd satisfies: 

Qd = (8.3) 

According to equations (|H3 and (IH3I), we have also: 

1 

- -g - 'i - 'd i-j—I 

We denote by 'I'i 'I'i {x,£) and 'I'± {x,£) the Bloch solutions described in section ETl We set: 

¥yix,q,,E)=^\x,E-W{q,)) ; E) E - W^q:)) ; 'i'Ux,q^, E) =W^\x, E - W{q})). 

We define the functions ujy, ujy and ujy associated by to the branches kg, ki and kd- 

8.1.3 Ideas of the method 

The computation is similar to this done in IHIEIEI It is based on some elementary principles that we 
outline now. 

1. Periodicity. 

The consistency condition lEHl) implies that the Wronskians are e-periodic in if. To get a total control 
of the Wronskians in a horizontal strip, we only need to control them in some vertical sub-strip of 
width e. 

2. Analyticity. 

Since the functions {<f,E) ^ f^{x,ip, E,e), {f, E) i—> ff{x,ip,E,e), {if,E) ^ fy{x,(f,E,e) are 
analytic on Sy x U, their Wronskians are analytic in {ip, E) £ Sy x U. This allows us to expand 
them into exponentially converging series. 

Let w(ip,E,e) be an analytic function in {ip,E) which is e-periodic in p. We set: 

w{f, E,e) = ^ Wk{E, £)e^ 
kez 

The Cauchy formula gives an estimate of the Fourier coefficients: 

Wk{E,£) = - / w{p,E,£)e ^ dp, Vfc G N, Vpo G Sy. (8.4) 

£ Jipo 

By moving Im pQ in [—Y,Y], we get a control of positive and negative coefficients. 

8.2 Asymptotic expansion of d{cp,E,e) 

In this section, we shall establish the following result. 

Proposition 9. For any Eq in J , there exist a complex neighborhood Uq of Eq and two functions 
{p, E, e) I—> b~ {p, E, e) and {p, E, e) b'^ {p, E, e) such that: 

• The coefficient d defined in llt).2(ill can be written: 

d{p,E,e) = i(J^w{fi,{fi)*)[b-{b+y - (6^)*6+]. (8.5) 

• The functions {p, E) ^ bg{p, E, e) and (jp, E) i—s- [p, E, e) are analytic on Sy x Uo- 
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• The functions <f bg {(p, E, e) and ip ((/j, E, e) are e-periodic and admit the following Fourier 
asymptotic expansion, when e ^ 0; 



bg{‘f,E,e) ='^{bg)k{E,e)e - , 

(8.6) 

with 

fcez 



{b~)Q{E,£) = aie~i-^o" (“+-“-)[! + o(l)], 

(8.7) 

and 

VfcfyO, l(bgY(E,e)l<Ce “/^e 'J 

(8.8) 


f’d(T,E,£) =J2(b+)k(E,£)e ^ , 

(8.9) 

with 

kez 

(b+)o(E,£) =iaiei-^°" + o(l)], 

(8.10) 


{b+YiE,e) =-ia^e^^r+ "fyf+ o(l)], 

(8.11) 

et 

VA:>1, \ib+Yi^,E,e)\<C\{b+YiE,e)\e 

(8.12) 


Vfc<0, \{b+)k{p’,E,£)\<Ce “/^e ' °, 

(8.13) 


The rest of the section is devoted to the proof of Proposition 
Fix Eq G J. According to the choice of Kg and k^ isections 15.9.21 and 15.9.It . there exist two analytic 
functions Oig{E) and ad{E) such that: 


(Og/l® )* = iagUghf, 

{adh%Y = iadadh'l. 

Now, we use the function fi constructed in Proposition |3 There exists a neighborhood Uq of Eq such 
that we can write: 

aghf = -icTg{bg)*fi + bgif,)*, 
adh% = -i(Jd{b'\)* fi + b'l{fiY. 

The coefficients ag and Ud are defined in equations and We compute: 


This leads to: 





ag{E)=e-iIo^ 
adiE) = , 


The coefficients bg and satisfy: 


bg = agUg , b+ = ada+, 
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where the coefhcients a and aj" are given by: 


and: 


We compute, for E gUq: 


aj = 


'w{fi,h^_) 

w{fzJt) ’ 

(8.14) 

'w{fi,h'l) 

HfzJf)' 

(8.15) 


d{ip, E, e) = w{agh^_,adh'l) 

= [K(^dT -^di^y] 

8.2.1 Continuation diagram of fi 

First, we describe the asymptotic behavior of the function fi in some domains of the complex plane. 

Lemma 8.1. We suppose that the assumptions of Proposition\^are satisfied. FixY < Y. Fix cpg < ipf 
and tpd > iff. There exists yo G]0, Im Lpi[ such that the function fi has the following asymptotic behavior: 
• For <p G {gj G Sy] Re if G YPg, Tr]}} fi ^o,s the standard asymptotic behavior: 


3 pjo 


<e-^^r‘^+ (^3 +o(l)) . 


fz = qge-^R 

For ip G {if G Sy', Re p G [pf, Pd]', I'm P > —ya}, fi has the standard asymptotic behavior: 

fi = iq^ef <e^rt (5-1 + o(l)) . 

For p G {p G Sy] Re p G [pf, Pd\] Im p < —yo}, fi has the standard asymptotic behavior: 


i IT S _ (271 fid} r‘^z 1 J _ ^ \ 

f. = -iqj^eBlo «-e ef° “+e (^'1 + o(l)) . 

Proof This lemma is similar to the continuation diagram presented in section 6 of |5]. Thus, we 
give only the main ideas of the study and refer to this paper for the details. The continuation diagram 
is represented in figure 1111 In this figure, the straight arrows indicate the use of continuation lemma 
iLemma 15.(111 . the circular arrows the use of the Stokes lemma ILemma 15.711 and the hatched zones the 
use of the Adjacent Canonical Domain Principle (Lemma^)). To complete the proof, it remains to explain 
the connections between the different objects of the WKB method. 

• According to the definitions given in section 15.9.21 the branches Ki and Kg are equal in S- and, for 
all p G S-, we have: 

Ki{p + 0) = Kg{p - 0), + 0) = TKt/j - 0), uj].{p + 0) = co^{p - 0) . (8.16) 

Besides, it remains to link qi and qg. Section [8.1.21 implies: 

VpGS-, qi{p)=qg{p). 

• Similarly, we have, for all p G S+: 

Ki(p-O) = -Kd(p + 0), 'i'±(x,p - 0) =’i'l^(x,p+ 0), (8.17) 

u;j.(p - 0) = u;l^(p + 0), qd(p + 0) = -i qi(p - 0) 

• We study finally the link between Ki and Kd along the Stokes line c beginning at We consider 

the quasi-momenta kd and ki associated to Kd and Ki. Equation for kd and fci, on either side 

of [E 2 , E3], implies that Kd and Ki satisfy the following relations, for p G c, 


. r-^ 

‘^+e (d^)! 


KdiF + 0 ) = 2 Tr - Ki{p - 0), ^'±(a:, p + 0 ) =p - 0) 

LO%{p + 0) = Lu'f{p - G) qd{p + Q) = iqi{p - Q) 


(8.18) 


❖ 


40 








Figure 11: Continuation Diagram 


8.2.2 Computation of and 

Now, we compute the coefficients b~ and 6^ given by and (Hna . 

According to Theorem 01 we know that the asymptotic behavior of the function h® remains valid in the 
domain {(p & Sy', Re (}p) G [tpg, Lemma IHtI gives the asymptotic behavior of ft in this domain and 

we get : 

yp&Sy, a~{p,E,e) =cri[l + o{l)]. (8.19) 

Fix Yq G]0,y[. In the strip Sy^, we write: 

ag{p,E,e) = '^ane^^ ( 8 . 20 ) 

The coefficients an satisfy: 

1 rVo+s 

an = - ttg {p, dp, Vn G N, V(/3o G {—Fq < Im < Fq}- (8-21) 

£ J(po 

Fix n > 0. We estimate |q;„|. We use formula (Unrii for Im po = —{Y — 6), and we get: 

lonl <Ce-2""(^-‘5)Ae^. 

We treat similarly the case n < 0 with Im po = {Y — S) and we obtain: 

\an\ < . 

Besides, we have: 

ao = cTi[l + o(l)]. 

We fix S For a constant C such that a < 27r(y — Fq) — i5(C + 27r), we obtain the estimates 

(18.711 and (18.8t . 

The arguments for the coefficients and are similar. 

8.3 Proof of Lemma 12.21 

Now, we want to express the coefficient d in a more understandable form. We begin with proving Lemma 
EH We recall that we denote by pf and pi, pi the branch points of the complex momentum, and by 
Er and Ei the related ends of tT(iFo)- We shall prove the lemma in the case (15.1511 . Let Ki be the branch 
described in section Em Ki satisfies (EHJ- 
We shall prove Lemma o for the branch Ki = Ki. 

• First, we express S and as integrals of the complex momentum along complex paths. Let 7 be 
an oriented curve, we call 7^ the curve oriented in the opposite direction. Fix G M and G M 
such that: 

Pd > pt ; Pg < Pr- 
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We define the complex paths 7$, 75 and 7^7: 


7$ = [if^ + iO, + iO] U — iO, — iO], 

7S = (o- + 0) U (a'^ - 0), 

7g,d = [^g + fO, 0 + iO] U {a+ - 0) U (cr^ + 0) U [0 + lO, ipd + *0]. 


These paths are represented in figure El We have the following result: 
Lemma 8.2. The coefficients <i>d and S can be written: 



/ 

J7. 



where k = Ki in S- and k is analytically continued along each path; k = Ki in S- and k is analytically 
continued along 7^. 


Proof 


— First, let us justify the fact that integrals along 7$ and 7s can be considered along closed curves. 
It suffices to show that k can be analytically continued along 7$ and 75. 

We consider the curve 7$. We have taken the cut of 7$ in (p~. We show that k has the same 
values on each side of the cut. k = Ki on [ipT + * 0 , ipf + iO], since k is continuous to the right of 
(pf, we obtain that k = —Ki on — i 0 ,(pf — iO]. In addition, K{(pf + iO) = 0 = K^pT ~ * 0 )i 
which proves that the integral can be taken on the closed curve 7$ . 

The arguments for 75 are similar. 

— We compute: 



/ 


1 

2 


Similarly, for the coefficient S{E) 



— It remains to study d)^. We introduce the branch m and we cut in elementary segments: 



= 2 



2^d{E) + 2Tr{ipd - Pg) 


This ends the proof of Lemma l ?01 0 
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Figure 12: Some complex paths 


• We use Lemma [8.21 to prove the analyticity of S and d>d- 

First, we consider d>. We can deform 7 $ to a closed curve going around and staying at a 

nonzero distance from this interval. Besides, n is analytic in E on the integration contour when E 
is close enough to J. The analysis of the coefficient S is done in the same way. To prove that is 
analytic, we deform the curves 'yg^d and 7 g_d to stay at a nonzero distance of the cross. 

• Fix E G J. On the interval the branch Hi satisfies Hi G [0, tt]. Thus, the function ^{E) is 

real positive on J. 

Now, we give a simplified expression of S: 


S{E) = -i 


/ {n — Ki{u))du + / {n — Ki{u))du = 2 Im / {n — Ki{u))du 

J ( 74 - ( 7 — J ( 74 - 


( 8 . 22 ) 


On cr, the branch Ki satisfies Ki G [ 0 , 7 r]. According to (12.811 and 118.2211 . we obtain that 0 < S(E) < 
27rlm ipi{E). 

Finally, we have {Ki{u) — Tr)du = (Ki(u) — 7r)du. Consequently, the coefficient is real. 

• Now, we compute S' and d)' on J. Let k be the branch of the Bloch momentum continuous through 
[Er,Ei], then K{ip) = k{E — W(}p)) and 


$'(£:)= f ' k'{E-W(u))du + k{E-Wiip+)) - k{E-Wiip-)) = j k'{E - W{u))du. 

J 


We recall that k has some branch points of square root type at the ends of spectral bands (see section 
consequently, the integral k'{E — W{u))du is convergent. In the interval [E^^Ei], k'{£) > 0 
and {Ei — Er)^' takes positive values on J. 

The analysis of S' is similar. 

• We complete this section with the following formulas: 


^d{E)+iS{E)= / Ki{u)du — 2 / —'!T){u)du + / Ki{u)du 

J J ( 7 _ J 

pQ p pO 

—^d{E) + iS{E) = — / Ki{u)du — 2 / (ni — Tr){u)du — / Ki{u)d 
J J ( 7 ^ J 

When = TT, the proof is analogous for the branch Ki = 2^ — Ki. 


(8.23) 

(8.24) 
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8.3.1 Further computations 

We recall that the functions and wL are defined in We consider the integrals of lo\ and uf‘_ 

along some paths of the complex plane. We have the following relations: 

Lemma 8.3. The integrals of uj\ and ujf satisfy: 


VEg j, 


E)du = 0, 


'[Vr ,vt] 


a;!_(u, E)du = 0 


JVpt ,vt\ 

'iE G j uj\{u, E)du = 0, I u!^_{u, E)du = 0 
J (7 J a 

There exists a real number p such that: 

VFeJ, [ {ujX{u,E)-ujUu,E))du- [ 

J[(pf, 0 ]Ucr^ J <7 

Proof We consider the case (EISJ. 


(8.25) 

(8.26) 


/cr_U[0,l/3r ] 


{pj\{u,E) — ijj'L{u,E))du = ip (8.27) 


We first prove According to E3, we compute: 

f E)du = — ( g\{E — W{u))W'{u)du = f 


E-W([Vr,V>^]) 


g\{e)de = 0 


Indeed, for E G J, the subset E — W{[ipf is a complex path of energies connecting E^ to 

Er and containing {E — IF(0)) g]Ei, E 2 [). We have shown this path in hgure IT!^. Particularly, 
E — W{[pf ,pf]) is a closed path and does not surround any pole of the meromorphic function g\. 
Consequently, the integral is zero. We prove similarly that 


j w!_(u, E)du = 0. 
"I .. 1 


• We consider now (18.2611. We write: 


f uj\{u^E)du = — f g\{e)de 
Ja JE-Wla) 


lE-W(a) 

The image of the path a is shown in hgure E)3. We deal with ujf similarly. 

• Finally, we compute: 

f E) — E))du — f {ijj\{u^E) — u>f{u^E))du 

J[Lpfg]Ua^ J{T-U[0,ipr] 

= [ {g+{e)-glie))de- f {g!^{e) - gl{e))de 

J E-W[[vt fi]yja+) d_E-W((T_u[0.<^y]) 

The images E — W{[pf , 0] U cr+) and E — W{cj- U [0, pf]) are two paths of energies connecting Er 
to Ei (see hgure EP)- By analyticity of {gf_ — gf) in the domain Re (e) G]Er, Ei[ , we obtain that: 

f {uj\{u,E) — uj'^_{u,E))du — f {ui\{u, E) — uj’_{u, E))du = 2 ( {g\ — g''_){e)de 

J [iff ,ei\VJ<7^ Ja-U[0,ipr] J Er 

(8.28) 

It remains to show that this coefficient is purely imaginary. To do that, we point out that (^h)* = g’f, 
according to KWi . Equation becomes: 


= 2 


/ — co^_){u, E)du — j {ijj\ — u}'‘_){u, E)du 

J<T^U[0,ipr] 

E ■ E' \ / E' E ■ \ 

^ g\{e)de-{gf.)*{e)de\ = 2 g!^{e)de - 5 +(e)dej 
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Er^ _ Ej 


Er 


A 


E, 


A 


B C 

Figure 13: Paths in the complex plane of energy 


This ends the proof of Lemma ESI ❖ 


8.4 Equation for the eigenvalues 

The following result gives a characterization of the eigenvalues of 
Proposition 10. We assume that {Hy), {Hw,r), {Hw,g) and {Hj) are satisfied. 

There exist Eq > 0, a neighborhood V = V of J, two functions {E,e) $(15, e) and {E,e) ^ $d(i5,£:) 

defined on Vx]0,£o[ and two functions {(p,E,e) i—*■ F{ip,E,e) and {ip,E,£) R 2 {lp^ E, e) defined on 
K X Vx]0,eo[ such that: 

1. E is an eigenvalue of if and only if: 


E{ip,E,£) = 0 


2. The function F satisfies: 


3. The 


e R, vt; e v, Ve e]o,eo[, 
function ip A((/3,15, e) is e-periodic and 


Its founer 


F{ip,E,e) = cos 


cos ( + ^+p\+ e-^(^)/-ll 2 ((^, 15, e) 


(8.29) 


4- The functions $, $d satisfy the following properties for any e g]0,£o[.' 

• 15 I— > $(15, £) and E i—> $d(15,£) are analytic on V. 

• $(15, e) = $(15) + o{e) and $d(15, e) = $d(15) + 0 (e) uniformly for E gV. 

5. For any e g]0,£o[, the function (ip,E) R 2 (ip, E, e) is analytic on K x V. Besides, there exists a 
constant a > 0 such that,for all e g]0,£o[, and all E in V, the function R 2 satisfies the following 
properties: 

/ l? 2 (n, 15, £)du = 0, / l?2(n,15,£)e = du = 0, / R 2 (u,E,£)e = du = 0, 

Jo Jo Jo 

sup |1?2(</5,15,£)| < e“¥ 

The functions $, $d, S are defined in T,em.m.n \2.‘A p is a real number defined in 

Now, we prove Pronositiou II01 
• Now, it suffices to compute the Fourier expansion of: 

bf(bd)*iT.E,£) =^jn(E,£)e^^. 
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By using the asymptotic expansion of the coefficients b and given in LemmalHl we prove that: 


7 o =-ie ’[l + o(l)]. 


H 7ne^ 
nez\{ 0 ,l} 


= 0(e pour (/j € • 


Actually, 


7 o = ao/3o + 2^ a«/3-« =-*e e L j[i + o(l)]. 

n^O 


According to 1123, we simplify: 



r^r 

/ {oj]. - oj'l)-\- 

/ 

Jo 

Jo 


= 0 . 


According to 123, r- = d’(A). Consequently, 


i^(E) 


7 q =-ie = [l + o(l)]. 

We compute: 

7l = aoPi + Unfil-n 
n^l 

We start with computing ao/3i. To do that, we deduce from equation that: 

y’O pO P 

/ Ki{(p)dip + / Kt{(p)d(p - / {Kt{(p) - Tr)d(p = ^diE)+iS{E). 

J J J (7 — 

Equation leads to: 

^0 ^0 p 

/ Ki{ip)d(p+ / Ki{ip)d(p- / {Ki{(p) - Tr)d(p = ^diE)+iS{E). 

J ip'^ J Pr ^ — 

Besides, according to Lemma IQ we have: 


rvt rO r 

/ u;\+ u;l- 

Jo J pr J CT-f 


— w!_) = ip- 


and: 


ao/3i = ie + o(l)]. 

Since S'(A) < 27rlm pi{E), we estimate the remainder in the expansion: 

I ^a„/3_„| = o{e~^>^). 

n^O 
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Finally, for p 0,1, we estimate: 


Tp ^ ^ ^nPp—n' 


For p > 1, we have: 


Similarly, we estimate for p < 0, 




o/ / 27 rYo(p-l) 

|7p| = -“-' 


C / / 2 TT Vq ( I P I — 1 ) 

|yp| = e-®/"e-“/"0(e--). 

• Now, we consider (p G K. We compute the Fourier asymptotic expansion of the coefficient d{E, ip, e) 
in a neighborhood Uo of Eq: 

^ 2zTi7r<^ —2z7i7rtp \ 

Xo{E,e) + ^ (A„(F;,£)e^^ + (A„)*(F;,e)e ^ ) j 

new / 

= i{wok’i){E - W (0)) E, e). 

neN 

where Un{p, E, e) = \n{E, E)e~^ + (A„)*(if, e)e , pour u G N*, et uo(p, E, e) = \o{E, £)■ 


We have: 


uo{(p,E,e) = 'yoiE,e) - i^iE^e) = -iE'^g{E,e) - ie * = p*(£',e). 


where g{E, e) = 1 + o(l). 

We define g{E, e) = rg{E, where the functions E rg{E, e) and E 0g{E, e) are analytic 

and satisfy 

^*g=^ 9 ^ ’'s = l + o(l) 0 *g = dg, 6»g = o(l). 


We simplify: 


Similarly, we compute: 


uo{ip, E, e) = -irg{E, e) cos ^ + (^g{E, e) ) • 


Ul 


{(fi, E, e) = irh{E, e)e cos + p + eh{E, e)^ . 


where the functions E ri^{E,e) and E 9fi{E,e) are analytic and satisfy 

rl = rh, rh = l + o{l) 01 = Bu, 9h = o{l). 

In addition, we have the following estimate of the remainder: 


Y^Up{p,Ep 

p>2 


-S(E) -a 

< Ce = e = pour G K. 


• We have proved that, for if in a neighborhood of Eq, the Fourier expansion of d{E,ip,e) can be 
written: 


+i|l + o(l)|e^ CO. (+ p + 0(1)1 + c-^O(c^). 


(8.30) 


The compactness of J implies that there exists a finite number of intervals {Jk\k^{i...p} such that: 
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1 . J c u Jk 

feG{l'p} 

2. For any fc G {1, • ’ ’ iP ~ 1}; the intervals Jk and Jk+i overlap. 

3. For any k G {1, • ’ ’ there exists a complex neighborhood Uk of Jk such that the expansion 

(non is satisfied on Uk ■ 

We shall prove that we can define some functions 3> and on the whole neighborhood V = 
IJ Uk- To do that, we shall “stick” the expansions obtained on each interval. 

feefi-p} 

The coefficient uq is written: 

VE G Jk, uo{E, e) = ro^k(E, e) cos + Qo,k{E, 

VF; G Jfc+i, uo{E,e) = ro,k+i{E,e)cos ^ + So,k+i{E, e)^ 

where ro,fc(F;,e) = l+o(l) and 0o,fe(£’, e) = o(l) (resp. ro,fc+i(^^,e) = l+o(l) and 0o,fe+i(Fl, e) = o(l)) 

for E G Jk (resp. E G Jfe+i).We get that: 

ro,fc(£;,£:) = ro,fc+i(£;,£:) = ro{E,e) et 6 o,k{E,e) = 6o,k+i{E,e) = Oo{E,e) for E G JkC] Jk+i 


The function $ defined by its restrictions to each Uk is analytic on V. 
The case of is treated similarly. 

Defining 

_ diip,E,e) _ 

i{wok'^{E-Wmro{E,ey 

we finish the proof of Proposition EH 


8.5 Localization of the eigenvalues 

In this section, we deduce Theorem ^ from PropositionllOl 
We solve equation E{ip, E, e) = 0, where F is described in (18.2911 . 


8.5.1 Energy levels E^’‘\e) 
For E G V, we start with solving: 


HE,e) 

cos-= 0 


(8.31) 


E I—!■ ^{E, e) is a real analytic function. For a sufficiently small eo) by Lemma 1^^ there exists a constant 
m > 0 such that: _ 

VEgV, V£:G]0,eo[, |$'(E,e)|>m (8.32) 

Consequently, equation IjOTll has a finite number of zeros in J. We denote them by E^ye)^ for I G 
{L-{e),L+{e)}. They are given by: 

$(F;W(e),e) 


= Itt ■ 


and satisfy: 


F;('+b(£) - F;(b(£) = 


VZG{L_(e),...,L+(£)}. 

1 


-7r£ + o(e). 


(8.33) 


(8.34) 


4>'(f;«(£)) 

The distances between two consecutive zeros are of order e. Precisely, by combining with (HOHl . 

we obtain that there exists a constant c > 0 such that: 


1 


-£ < \Ed+^\e) - E^^\e)\ < ce, Wl G {L_{e),..., L+ie) - 1} 


(8.35) 


First, we prove that the zeros of F are in an exponentially small neighborhood of the points E^ye). 
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8.5.2 First order approximation 

We give a first order approximation of the zeros of F. 
We set 


ao{E, e) = cos 


HE,e) 


We can assume that the neighborhood V is sufficiently small and such that, for any E gV, 


Re {S{E)) > (3>0. 


Then, there exists a positive constant A such that 

\F{(p,E,e) - ao{E,e)\ < 
In addition, we have the following inequality: 


3C > 0/ 


^E,e) 

cos- 

e 


>jd{E, U £;«(£)). 


Actually, there exists a constant c > 0 such that: 

I cos 9\ > cd{9, ttZ + tt/2 ). 

By using we obtain the relation: 

|$(A,e)-$(A(')(e),e)| > m\E - E^^'> {e)\ 


and finally: 


cos • 


HE,e) 


>|dlA, U A«(£)|. 

' -.L+} 


For zo G C and r > 0, we define 


(8.36) 


D{zo,r) = {z gC ■, \z-zo\<r}. 

Inequality 18.3611 implies that there are no zeros of F outside exponentially small neighborhoods of the 
points E'd3{£). Precisely, there exists a positive constant D such that, if r > then for any 

E G dD{E'd){£)^r), we have: 

\F{ip,E,e) - ao(A,e)| < |ao(A,e)|. 

Rouche’s Theorem implies that, for any I, E has exactly one zero Ei{(p,£), in each neighborhood 
D{Ed'> [e), Dee~^^'^) of The relation E = F* allows us to recover that the eigenvalues are 

real. Indeed, if F{E) = 0, E is also a zero of E. By uniqueness, we obtain that E = E. 

We set: 

Ei{(p,e) = Ed\e) +£\i{ip,e). 

We know that \i[tf,£) is exponentially small. Now, we compute its asymptotic behavior. 

8.5.3 Second order approximation 

We define: 

ai{(p, E, e) = F{lp, E, e) - ao{E, e). 

We write 

g-S(£,(v,E))/e ^ ^ 0(Ai(^,e))). 
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Similarly, with the help of the modified phase we obtain the expansion: 


cos 


^d{Ei{(p,e)) + 27Tip + pe \ ^ f ^d{E^^\e)) + 2Tr(p + pe \ 


The expansion of oi can be written: 

ai{p, Ei{(p, e), e) = ai{(p, (e), e)(l + r((/7, (e), e)). 

Moreover, we use the first order Taylor’s expansion of the function E ao{E,e): 

ao{Ei{^, s),e) = (-l)'+i$'(if«(£), e)(l + r{^, , e)) = (-l)'+i<i>'(ifW)AK(^, e)(l + o(l)). 

By combining these computations, we finally obtain: 




(- 1 ) 


i+i 




^-S{E(‘He))/e 


COS 


^diE^'‘\£),£) + 27r(^ 


+ P +o(l) • 


8.6 Application to the trace formula 

In 0, the author proves the existence of an asymptotic expansion of tr[f{Hy,^e)], for / G , when 
Supp / is disjoint from the bands of E[q] in addition, he computes explicitly the first and second terms 
of this expansion. 

Corollary fallows us to recover these terms. 


8 . 6.1 


Let J be an interval satisfying (iLj). Particularly, J is such that J C (ctqc U a sc) = 0- For / G C^, with 
Supp / C J, we compute: 

tr[/(iL^.,)] = ^ fmv,£))- 

l£{L- 

Let /3 > 0 be such that S{E) > (3 for any E G J] according to Theorem Q] we know that there exists a 
constant C > 0 such that: 

VuG[0,£], \tr[f(H^^c)]-tr[f{Hu,e)]\ <C ^ 

iGfi- (E),...,r+(e)} 

By integrating with respect to u, we obtain that: 

ti- [fiH^,e)] = -[ tr [f{Hu,e)]du + 0{e~^^^). 

£ Jo 

According to Theorem we know that there exists a constant C such that 


Vm G [0,£] 


tr[/(iJ„,,)] - ^ /(A«(^)) 

/e{L-(e),...,L+(e)} 


< 




By integration, we obtain: 

1 


r tr [f{Hu,e)]du = Y. (^)) + 0(e-^/^) 

^0 ,^fr /_X r /-M 




Now, we estimate: 


^ f{Ed){e))= Y fo^-\e{ln + 7 T/ 2 )) 

iGfi- (e),....r+(e)} lG{L-{e),...,L+(e)} 


(8.37) 
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8.6.2 

Now, we compute this last term. 

Lemma 8.4. Let f be a funetion in (7“ sueh that Supp f G J. The trace of has the following 
asymptotic behavior: 

r tr[f{H^^,)]du = - [ f{£W{£, e)d£ + 0{e^) 

Jo Jj 

Proof The proof of this Lemma is based on elementary results of real analysis. 

• We apply the Poisson formula to the function / o e 

sJ2fo ^-\e{l7r + 7r/2)) = 2 ^(-l)-(/T^i) . 

iez raSZ V ^ / 

Besides, the Fourier transform of / o satisfies the estimates: 


Vi/ > 1, 3 Cj/ > 0, such that 


(/vr.) (^) 


e 

< a—. 


Actually, since / o $ ^ is f o is bounded. 

This leads to: 

£ 'Y2 f ° + 7r/2)) = 2(/ o $“i)(0) + 0{e°°). 


• It remains to prove that: 


2(/o$-i)(0) = - //o$-i(M)du. 


With the substitution u = 4‘(if’), we obtain that: 


2(/o<l>-i)(0) = - / fi£)^'{£,£)d£ 


This completes the proof of Lemma EH 


❖ 


8.6.3 Conclusion 

To get an asymptotic expansion of the trace at any order, it suffices to know an asymptotic expansion 
of the modified phase at any order. Our computations are not accurate enough, but we know that 
<!)'(£) = + o{e), hence: 

i / fi£)^'{£)d£ = - f fi£)^'i£)d£ + oie). 


To transform the right member of previous equality, we do the substitution (k, u) i—> {E{k) + VF(m),u), 
which implies: 


- [ f{£)<i>'i£)d£ = ^ [ r f{E{K) + W{u))d^du 

J J t/[—TT jTt] j ifr 


We finally obtain: 


[ tr [f{Hu,e)]du = ^ [ [ f{E{K)+ W{u))dKdu + o{£) 

Jo i/[— tTjTt] j ip~ 


This ends the proof of Corollary Q] 
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8.7 Asymptotic behavior of the eigenvalues 

Now, we give a second application of Theorem ^ for the computation of the asymptotic behavior of the 
eigenvalues of Such a computation is outlined in |3], in the case V = 0. We obtain an explicit 

result at first order. 

Under the assumptions of Theorem^ is the only end of (7{Hq) belonging to {E — 1U)(]R). We define: 


d>p{Eji^ 


lim 


E^Bn, 

[-E'n 


k{E) - fc(U^) 
\/E — En 


(8.38) 


Corollary 2. Let verify the assumptions of Theorem^ The eigenvalues E^’'\(p,e) of Elhave 
the following asymptotic behavior: 

E^‘\g,,e) = ^-\e{l7r + 7r/2)) + 0{e^). 

Particularly, E^^'> [ip, e) has the following Taylor expansion at first order in e : 


E<^^\p,e)=Er + W{0) + 



1 

di[Er) 


[21 + l)e + o[£), 


where di[Er) is defined by (18.3811 . 

Proof The first equality is obvious. It suffices to give an expansion of $”^(e:(/7r + tt/2)). To do 
that, we compute an expansion at first order of: 

(£'r + VF_+a) 

^{Er + W{0) + a) = / k{Er + W{0) + a - W{u))du. 

Jifr {Er-\-W- +Ck) 


The mapping IT is a bijection from [0,v5i^] to \W[Q),Er\. By the substitution av 
we get that: 


nvf(Er+W(0)+a) 


■f 


k[Er + av) 


W oW-^[W[Q)+a[l-v)) 


dv. 


T3,,+ 1™ k{E,.+av) _ 

a™ WoW-^{W-+a(l-v)) 

Similarly, on [</?“, 0], we have: 


di{Er-) y/v 
^2lV”(0) \/l-« 


W[Q) + a-W[u), 


^[Eg + 1U(0) + a) — di[Er) — ^ 

Consequently, by inverting the expansion of <1> in the neighborhood of Eg + W(0), we prove the result. 0 
We point out that, as in Ei a more accurate asymptotic expansion of $ would give a better result on 
the eigenvalues. 
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